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1 Introduction 


Quantum technologies promise to revolutionize the world we live in by using 
devices that utilize quantum effects [1] such as quantum state supersposition 
and entanglement to provide access to features that go beyond what is classically 
possible. With the 2012 Nobel prize in physics awarded to David Wineland and 
Serge Haroche [2] "for ground-breaking experimental methods that enable measuring and 
manipulation of individual quantum systems”, the expectation on quantum technology 
is now higher than ever. For example, the impact on quantum chemistry by 
simulating complex molecules is probably one of the most anticipated applications 
of quantum simulation [3], while encryption software companies are already 
developing and implementing post-quantum cryptography algorithms to counteract 
the computational speedup [4]. 


Billions of Euros of public funding worldwide, such as the multi-national European 
quantum technology flagship initiative [5], illustrate not only the technological 
importance of quantum technologies for the future, but also the hope for a long-term 
economical impact similar to the mid-20th century semiconductor industry, which 
resulted in the success story of the Silicon Valley. At the same time, large international 
companies such as Google [6], IBM [7, 8], and Microsoft [9] invest heavily in their 
own quantum research groups, attracting the attention of scientists, the media, and 
the general public alike. 


In analogy to classical computers, quantum information processing builds on 
quantum bits (qubits) as its basic processing units [10]. While a classical bit operates 
in binary 0 and 1, a qubit operates in a quantum superposition state |Y) = a |g) + B |e) 
that when measured collapses into either its ground state |g) or excited state |e). 
Although a single qubit does not offer a particular advantage over a classical 
bit, combining multiple interacting qubits in a joint Hilbert space potentionally 
offers exponential computational speedup compared to a classical computer for a 
certain class of problems [11]. A recently published paper by the Google group [12] 
claims to have shown the proposed quantum supremacy for the first time ever by 
performing a computational task on a 53 individually controllable qubit chip, which 
"for a state-of-the-art classical supercomputer would take approximately 10,000 years.”, 
with the significance of this archievement compared to the Wright’s brothers first 
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flight: “It is what the event represented, rather than what it practically accomplished, that 
was paramount” [13]. 


Of the different technologies nowadays considered as the basis for quantum 
computing, qubits based on superconductivity [14] are currently at the forefront 
of a solid state approach [15]. Superconducting nanostructures that utilize the 
Josephson effect result in a non-linear and non-dissipative circuit element, which 
forms the base building block of such a qubit. Together with quantum versions 
of capacitors and inductors, quantum electrodynamic architectures combine these 
qubits with superconducting resonators to encoded, manipulate, and retrieve 
quantum information using microwave photons [16, 17]. 


Despite the engineering progress and technological developments in recent years, 
there are still certain limitations in superconducting quantum circuits. Typically 
operating in the single digit gigahertz regime, such a system requires operation 
in a cryogenic environment at temperatures near absolute zero to freeze out 
thermal excitations and minimal probing power to not suppress the quantum 
effects. Qubit parameter fluctuation [18-20], unintentional crosstalk between on-chip 
components, and the need for constant elaborate error-correction while performing 
algorithms due to insufficient coherence times [21] renders the implementation of a 
universal superconducting quantum computer quite challenging. On top comes the 
connection to the outside world and the need for efficient transport and storage of 
quantum information. 


Intrinsic advantages and disadvantages of all technological quantum platforms 
lead to the emergence of quantum hybrid devices, which ultimatly combine the 
advantages of the nature of each different physical system [22] and aim for quantum 
communication and quantum enhanced sensing. Nitrogen vacancy based spin 
systems [23] with their large coherence times, for example, are hybridized with 
photons to convert quantum excitation into the optical photon band band for 
quantum network communication [24, 25] and nuclear spins for quantum memories 
exceeding over one second [26]. 


The engineered level structure of a superconducting qubit, basically an artificial 
atom with designed energy transitions, makes it ideal for the study of the interaction 
of matter with other physical systems beyond natural limitations [27] and quantum 
sensing applications [28, 29], which rely on the high sensitivity of any coherent 
quantum system regarding external influences to precisely measure physical 
quantities. The sensing of small magnetic fields [30, 31] or the coupling to single 
electron or nuclear spins with superconducting qubits is of great interest for quantum 
application. On the other hand, the direct coupling between the spin of a single 
electron and a superconducting resonator is too weak to reach the coherent coupling 
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regime. It is therefore necessary to enhance the coupling strength by using a large 
number of spins in paramagnetic spin ensembles [32], i.e. based on erbium doped 
crystals [33-35]. Recently, the scope of quantum hybrid system was expanded to 
collective bulk magnetic excitation. 


The general phenomenon of magnetism has been known for a long time in human 
history and is omnipresent in modern everyday life beyond the field of science. 
The physical origins of bulk magnetism in solids however, are purely quantum 
mechanical as stated in the Bohr-van Leeuwen theorem [36, 37]. The exchange 
interaction between magnetic moments in a ferromagentic material allows for 
collective spin excitation, phase-coherent precession of the magnetization vector [38, 
39], that spread over the crystal lattice and are described as a wave similar to 
phonon excitation. These spin waves, quantized as magnons [40, 41], are utilized 
in ferromagnetic resonance (FMR) experiments, which are used as a powerful 
technique for investigating magnetic properties in materials. The excitation of spin 
waves by microwave photons via magnetic dipole interaction [42] or Brillouin Light 
scattering [43], as well as the macroscopic quantum phenomenon of Bose-Einstein 
condensation of magnons [44] stimulated the field of magnonics [45, 46], aiming at 
transport and information processing via magnons on the nanoscale. 


The concept of a strongly coupled spin-photon hybrid excitation using ferromagnets 
was proposed [47-49] and realized [50] using gallium doped yttrium iron garnet (YIG), 
with the hybridized modes forming a cavity magnon polariton quasiparticle [51, 52]. 
With their large spin-spin coupling and high spin density, ferromagnetic excitations 
are described by a macro spin, which results in a higher magnon-photon coupling 
strength compared to paramagnetic spin ensembles. Subsequent cavity magnon 
polariton experiments using three dimensional microwave cavities explore the 
system at room temperature [53] and show ultra strong collective coupling [53-55], 
with the general concept also being expanded to cavity optomagnonics using optical 
photons [56-59] or mechanical degrees of freedom [57, 60]. 


Strong coupling approaching the quantum limit of excitation [61-64] opens the door 
to nonlinear hybrid quantum systems containing magnons and superconducting 
qubits [65, 66], using the microwave photon as common coupling partner, see 
Fig.1.1. With all the previous systems being purely harmonic, this allows for 
the first time to address the quantum properties of single magnons and opens 
the field of quantum magnonics with strong and coherent interaction between 
superconducting qubits and magnetostatic modes. With further studies counting 
magnon excitation [67], single magnon detection [68], and the quantum simulation 
of quasi-particles consisting of fermions and bosons [69], magnon-based hybrid 
quantum system offer a large opportunity for new quantum technologies such as 
quantum sensing [61, 70], quantum transducer [71] or the search for dark matter via 
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axion to magnon conversion [72-74]. The main challenge hereby is to mediate the 
opposing nature of magnetic fields and superconductivity, with the need for spatial 
separation and a quantum bus in a three dimensional approach, or magnetic fields 
confined in a small volume in two dimensions. 


conducting 
Qubits 


Figure 1.1: Quantum magnonic coupling scheme Quantum magnonics combines collective magne- 
tostatic modes with superconducting qubits via indirect photon-mediated coupling using different 
components of the electromagnetic field. Inspired by [75]. 


In this thesis, we study nonlinear quantum circuits based on magnetic materials, with 
the focus on yttrium iron garnet due to its unmatched low internal damping and 
high spin density. As the basis for this thesis, we introduce collective spin excitation 
as a harmonic system and describe the spin-photon coupling mechanism with 
the focus on cavity-magnonics, where microwave photons in a three dimensional 
cavity resonator couple with magnons via their magnetic field component. We then 
establish the superconducting quantum bit as a nonlinear element and discuss how 
a qubit is operated in a quantum circuit. The following chapter combines the two 
systems under the concept of quantum magnonics, where we introduce effective 
coupling mechanisms between magnetic excitation and the qubit via microwave 
modes. A chapter about the experimental methods used in this thesis and basic 
resonator characterization measurement prepare for the measurements of various 
hybrid systems. 


The experimental part of the thesis is divided into two main chapters. First, a 
strongly coupled harmonic magnon-photon system is studied at the experimental 
working conditions of a superconducting quantum circuit. We probe the system 
using continuous microwave spectroscopy, with the focus on the magnetic resonance 
linewidth in order to identify the dominating loss mechanisms of magnetic excitations 
at various experimental conditions from several kelvin [76] down to millikelvin 
temperatures and single excitation powers [62]. In addition, we focus on the influence 
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of the energy mixing ratio between the magnon and photon share in the resulting 
cavity magnon polariton quasiparticle on the magnetic linewidth. 


In the second experimental chapter we add the superconducting qubit to the system, 
containing now three circuit elements. The nonlinear quantum system is studied in 
two opposite regimes, which mainly differ in the frequency detuning between the 
qubit and the microwave frequency. In the dispersive regime, at a large detuning 
on the order of 100 megahertz, we study the coherent interaction between single 
magnetic excitation and the superconducting qubit, where we demonstrate indirect 
qubit-magnon coupling and coupling enhancement by the participating number 
of spins. We modify the cavity resonance to push the system to the opposite 
end of the detuning spectrum, where qubit and microwave photon are almost on 
resonance, with a detuning on the order of ten megahertz. We spectroscopically 
study the coupling strengths of the tripartite system and find coherent coupling 
at two working points. Thermally excited multi-photon transitions are observed 
at various temperatures and identified by simulation. The thesis concludes with 
time dependent measurements regarding Rabi oscillations on the tripartite system 
as a prerequisite for more sophisticated implementations of magnetic excitation in 
nonlinear quantum circuits. 


2 Cavity electromagnonics 


This chapter presents the basic building blocks of interaction between magnetic 
excitations (magnons) and microwave cavity photons, which is mediated by the 
magnetic field component of the photon being the connecting component. We 
approach the hybrid system from the side of cavity electromagnonics, where photons 
at a fixed resonator frequency with confined magnetic fields are used to probe 
magnetic excitations. At matching frequencies, the coupled system is described by a 
quasiparticle, the cavity magnon polariton. This stands in contrast to traditional 
ferromagnetic resonance (FMR) experiments, which are usually performed with 
photons in transmission lines at different frequencies, ranging over several GHz. 
Despite being able to reduce the excitation number to single digits by reducing 
thermal excitation, there are no quantum effects present in the coupled system, 
which in our case is purely harmonic. We begin by discussing microwave resonators 
and afterwards introduce collective spin excitations in magnetic materials which are 
described as quasiparticles that represent the bosonic states of a harmonic oscillator. 
We continue with the coupling between the two systems and their coupling strength, 
as well as the microwave reflection spectrum of the coupled system. The chapter 
concludes with a discussion about the used magnetic material and different magnetic 
loss mechanisms. 


2.1 Microwave cavity resonators 


We probe the magnetic excitation with microwaves and use the magnetic field 
component of a three dimensional microwave resonator to excite the magnons. 
The amount of excitation energy expressed in number of photons is an important 
experimental parameter, since many effects in the quantum regime are power 
dependent. We estimate the average photon occupation in the resonator by its 
characteristic quality factors that describe its ability to couple excitation in and out 
and pass it to the magnetic sample. 


A three dimensional microwave cavity resonator is formed by a closed box, where 
modes of standing waves of the electromagnetic field form the resonance photon. The 
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energy is stored within the cavity and dissipated in the metallic walls or the filling 
dielectric material. While a cavity can hold an infinite number of resonance modes, 
there is a dominant mode given by the lowest resonance frequency. Resonances in a 
rectangular cavity are described by their wave number 


mm y? nT \ 2 Im\ ? 

Kmni i 4 ) j ( b ) j ( d ) , (2.1) 
with m, n, and | describe the number of anti-nodes in the standing wave pattern and 
a, b, and d the edge length in the x, y, and z direction, respectively. See Fig. 2.1 for a 
depiction of a rectangular cavity and the electric field components of the transverse 


electrical (TE) modes TE¡p1 and TEj92. Our notation here is the same as in Ref. [77]. 
The resonance frequency of all modes is given by 


A (EY (NPC ay 


with c being the speed of light and yu, and e, the filling material's properties. 


a b 


0.04 


0.02 


0.00 


—0.02 


3m[S;1] (arb. units) 


=0.04 


-0.06  -0.04 -0.02 0.00 
Re[ Si] (arb. units) 


Figure 2.1: Rectangular cavity resonator (a) Electric field variation. The electric field is parallel to the 
shortest edge of the rectangle, with the number of anti-nodes in each direction indicating the mode. In 
this case b < a < d, the dominant mode is the TEıoı mode. The magnetic field components circle in a 
perpendicular plane around the electric field maxima (not shown). (b) Cavity resonance spectrum and 
corresponding fit to the scattering parameter Sj; in the complex plane measured in the single photon 
limit with —140 dBm power applied at the SMA port. From the fit we obtain a resonance frequency 
w,/2rı = 5.239 GHz, and Q factors Qı = 3084.4, Q. = 5438.9, and Q; = 7124.9. The resonator is used 
for the experiments in Sec. 7.2. 


The cavity resonances are excited either capacitively to the electric field using the 
inner conductor of a coaxial cable, or inductively by a coupling loop penetrated 
by the magnetic field component. For resonances in the GHz regime the cavity 
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dimensions are in the range of mm to cm, large enough to store for example a 
magnetic bulk sample inside. Fabricating (and accessing) the cavity requires a cut in 
the design. This cut is made in parallel to the electric field to allow an undisturbed 
flow of electrons in the metallic walls. 


2.1.1 Quality factors 


The quality factors (Q factors) describe the losses of a resonator at the resonance 
frequency wo. In general a quality factor is defined as the ratio between the stored 
energy within the resonator Etot and the energy loss rate Pjoss, 


Etot _ wo 
Ploss 2K 


Q = wo 


(2.3) 


Here, « is the half width at half maximum (HWHM) of the resonator’s frequency 
response. There are three different quality factors in a resonator: The total loss of the 
resonator is described by the loaded quality factor Q}. It is defined as the reciprocal 
sum of the coupling quality factor Q. and the internal quality factor Qi, 


QT! = Qt Oe: (2.4) 
With Eq. (2.3) we find for the corresponding linewidths, 


K] = Ke + Ki. (2.5) 


2.1.2 Scattering parameters 


We describe the measurement signal of our resonator by a scattering matrix that 
contains the reflection coefficient r (w) of a probe signal with frequency w to our 
resonator, which is capacitively coupled via a coaxial cable to the outside world. The 
components of the scattering matrix for any generic multiport microwave network 
relate the amplitude and phase of the signal reflected from port i to the incoming 
signal of port j, 
yout 
Si pm ar . 
J 


(2.6) 


Our resonator has only one port, the reflection coefficient here is the ratio between 
the outgoing and the incoming voltage signal, 


Vout (w ) 


r (w) = Sı = V; (w) . 


(2.7) 
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We model a lossless resonator by an impedance Z, that is coupled to an open-ended 
transmission line with a characteristic impedance Zo describing the coaxial cable. 
Following Refs. [77-80], the charge of the input capacitor plate of the probed 
resonator is described by a harmonic oscillator that is driven by the incoming waves 
and damped by the transmission line with damping rate 2x. = Zo/L;,, with Ly 
being the resonator inductance. We find the reflection parameter depending on 
both impedances 

= Zr (w) = Zo 
a Zr (w) +Zo i 


The resonator impedance near the resonance frequency wr is approximated to be 


r (w) (2.8) 


Zr (w) = 2jLo (w — wr) (2.9) 


and we find the reflection coefficient 


w — Wy +jKe 
w — Wy — jKe 
2Ke 
j (w — wr) + Ko 


r (w) = Su (w) = (2.10) 


=1 (2.11) 


The internal losses in the cavity x; are accounted for by a complex renormalization 
of the resonator frequency wr —> wr + jx; [81], leaving 


2Ke 2Ke 
: =1-; 
j (w — wr) + Ki + Ke j(w= wr) +K 


$1 =1 (2.12) 
This allows us to fit the experimental resonance data and characterize a resonator 
by its resonance frequency and Q factors. We fit the complex reflection coefficient 
with a circle [82, 83], see Fig. 2.1 (b). The derivation based on the circuit impedances 
employs the usual electrical engineering sign conversions with the imaginary unit 
j = —i [79], in contrast to the derivation by the input-output formalism which 
uses the physics notation with the imaginary unit i. We will use the input-output 
formalism later to describe our coupled systems (Secs. 2.3.1 and 4.2). 


2.1.3 Average photon population 


The Q factors allow us to estimate the average photon population (n) in the resonator 
by a general power balance of incoming, outgoing, and lost power, 


Poat =: Pin Poss (2.13) 
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With the absolute squared scattering parameter defined to be |S11 |2 = e and the 


internal Q factor Qi = wo pa = hws ee we find 
= Ploss = l= ISu]? 
(n) = Qi ho? QiPin ie (2.14) 
Evaluating |S, E on resonance we find 
QF Pin 
(n) = Oc fiw? (2.15) 


Pin here is in units of watt, in contrast to the usual dBm used throughout the thesis. 
Converting the power in Eq. (2.15) into the power level Lp = 10log;p (rw) with 
units dBm we find 


2 
(n (Le) = 20 -10-6W. (2.16) 


2.2 Collective magnetic excitations 


Magnetic excitations in macroscopic ferromagnetic sample spread across the whole 
sample and manifest as collective excitation. Similar to phonons they are described 
as waves. Instead of exciting the crystal lattice, magnetic excitations form spin 
waves that describe the dynamic eigenmodes of the magnetization in a magnetically 
ordered medium [84]. Classically, spin waves represent collective, phase-coherent 
precessional motion of magnetic moments, see Fig. 2.2 that are quantized with their 
associated quasiparticles called magnons [40, 41]. In the magnetic ground state 
all spins in the lattice are in parallel. While an excitation can be assumed as one 
spin at one lattice site j being flipped, this is not an eigenstate of the Heisenberg 
Hamiltonian for quantum mechanical magnetic exchange interaction, which is based 
on the Pauli principle. Instead of a single spin flip, a lower-energy delocalized state 
is excited that is described by a linear combination of all possible spin flips. We 
describe an ensemble of N; spins in a static external field B || ez by 


Ns 
A =-—g BB, Y 5 -2J Y, $;: Sj, (2.17) 
i (i,j) 


where g* is the g-factor, ug the Bohr magneton, and $; the spin operator of the 
ith spin [85, 86]. The first term describes a dipolar Zeeman coupling and the 
second term the nearest neighbor coupling with exchange coupling J > 0 in the 
ferromagnetic case. 
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We employ the Holstein-Primakoff transformation [40] to map the spin operators Si 
to bosonic creation and annihilation operators m! and m, quantifying the number of 
magnetic excitation quasiparticles. Magnons are bosonic eigenstates of a magnetic 
quantum harmonic oscillator with the number of magnons corresponding to the 
reduction of the z-component of the total spin (see Fig. 2.3), resulting in a maximum 
of 2s magnons. We find the spin component operators 


St = SY +18) = 1/28 — fmit fu mj, (2.18) 
> = 87 18) = mt /2s — mth, and (2.19) 


Im, (2.20) 


with spin s on each lattice site i and commutators IM; m] = 6; [87]. 


AITTTETTTTT 


Ka 


b 


COO OOOO OTS 


Figure 2.2: Spin wave on a line of spins The spin excitation spreads over the crystal and all spins 
precess around their equilibrium position forming a spin wave. The picture shows an exemplary spin 
wave with one wavelength in perspective (a) and from above (b). 


Inserting the spin operators into the ferromagnetic Hamiltonian, taking all terms 
into the quadratic order of the bosonic operators, and Fourier transforming them 
into momentum space simplifies the ferromagnetic Hamiltonian to 


Ê = Y ham (k) iy, (2.21) 
k 


with the wave vector of the spin-wave mode k and the magnon dispersion relation 
Wm (k) [85, 86] with neglected constant ground state energy Ey. Considering only 
nearest neighbour interactions and a simple cubic lattice simplifies the dispersion 
relation in the long-wavelength limit to [66] 


ħwm (k) © g*ugBz + 2sJaż |k|?, (2.22) 
with the lattice constant ag. The contribution of the exchange interaction becomes 


negligible in the limit ag |k| << 1 and the magnon modes correspond to magnetostatic 
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modes generated by magnetization oscillation obeying the Landau-Lifshitz-Gilbert 
equation 


A ay w), (2.23) 


Fr Mx Hett = Me x Ar 
with the magnetization M, the saturation magenetization Ms, the effective field Hey, 
the electron gyromagnetic ratio y, and the Gilbert damping factor «. For spheroidal 
samples the modes were solved by Walker [88-90] and in the simplest Walker mode, 
the Kittel mode [91], all spins precess uniformly in phase with the same amplitude. 


This reduces Eq. (2.21) to 


A = ħwm ñ, (2.24) 
* 
with wm = È = Ba. (2.25) 
zZ 
(mim) =0 (S2) =s 
(mtm) =1 (Sy =s=1 
(mim) =2 ($2) =5-2 
(mim) = 2s (S2) = -s 
x y 


Figure 2.3: Holstein-Primakoff transformation Using the Holstein-Primakoff transformation the 
collective magnetic excitations are described as bosonic quasiparticles, which are distributed over the 
whole lattice. The number of magnons (mitm) corresponds to the reduced z-component of the spin 
quantum number. Inspired by [66]. 


A Kittel mode magnon corresponds to one single excited spin 1/2 particle in the 
field Bz, magnons in the macroscopic sample are therefore approximated by a 
single macrospin. The quantum harmonic oscillator approximation of Kittel mode 
magnons in the magnetostatic limit is only valid if the number of excited magnons 
is much smaller than the total number of participating spins, comparable to a 
small-angle approximation. We are able to tune the magnon frequency by applying 
a homogeneous magnetic field wm « Bz, and excite the magnon by a magnetic 
field oscillating with wm perpendicular to Bz. A non-uniform ac field results in 
higher order magnetostatic modes being excitable as well, which can be excited eiter 
voluntarily [54, 70] or as an involuntary effect [53, 61, 62, 65, 76] together with the 
Kittel mode. 
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2.3 Magnetic dipole interaction: cavity and magnon 


We describe the interplay between a magnetostatic mode and a microwave cavity 
mode by a magnetic dipole interaction. The cavity resonance modes are expressed 
by harmonic oscillator modes with resonance frequency w, and creator /amnihilator 
operators ât /ä, 

Ae = hw,a*a. (2.26) 
The interaction with the magnetic spins is given by a Zeeman-type 
interaction Hamiltonian 


Ns 
m-ce = 8* HB )_B (ri) $; 
i 


AN 1.8 (apat 
= 313) OB (ri) - $; (a+a Y (2.27) 
i 


with the microwave cavity magnetic field at position r; of spin i, B(r;) = 
ôB (r;) (a+. Using the Holstein-Primakoff transformation, we replace the 
spin operator by the sum of of magnon operators í1* and ñ multiplied by their 
orthonormal spatial mode functions sy (r;) [66], 


$; = V2s Y Sn (r;) L > : (2.28) 
n 


With the coupling strength g much smaller than the resonance frequencies we apply 
a rotating wave approximation and find 


Amc = eya F I dr 5B (r) -sn (r) (am! + atm) (2.29) 
= yng) (art + Atr) ; (2.30) 


where V is the sample volume and g”) the effective coupling strength between the 
magnetic mode n and the cavity mode. 


In the special case of a microwave magnetic field ôB (r) uniformly distributed over 
the magnetic sample, only the Kittel mode leaves a non-vanishing contribution to 
the integral over the N; spins within the sample volume. Using the average single 
spin coupling strength ¢ for all participating spins we find a YN, enhancement 
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compared to the single spin coupling, $ = \/Y; |g (ri)|? / Ns [92, 93]. The interaction 


Hamiltonian and the coupling strength are given by 
Am-c = hg (am + atm) and (2.31) 
c= sub a [SB (r)| V2Nss, (2.32) 


with the amplitude of the vacuum fluctuations |B (r;)|. This Hamiltonian represents 
the Tavis-Cummings model, which describes the interaction of multiple identical 
two-state systems with a single quantum field mode [94]. The collectively enhanced 
coupling strength accounts for the fact that a single magnon excitation is shared 
between all participating spins. The coherent superposition of the N, identical 
spins are called Dicke states [95]. The magnetic field at point r in the cavity mode 
is given by [96, 97] 


B.=i|8B (r)| (f(r)a— ft (1) 2), (2.33) 
with a dimensionless vector function f describing the spatial structure of the 
magnetic field mode. This gives a cavity ground state energy 


(0| ic |2 ro) = nr (2.34) 
hw 
(0| = 18B[? Vnoael0) = A (2.35) 
Ho 
SB] = , | Hote (2.36) 
2 V mode 


with the effective magnetic mode volume Vmode = f |f (r)|? Pr. The overlap 
between the magnetostatic and the electromagnetic microwave mode is accounted 


for by an additional factor that y = y + B2) /max { 8?) [53, 98], with the xz 
plane being in parallel to the cavity cut described in Sec. 5.3. The magnon-photon 


coupling is given by 


= den, POE ange. 2.37 
2 Vox = 


with the gyromagnetic ratio of the electron Ye = g* p/h. 


Similar to the pure cavity resonator measurements we probe our coupled magnon- 
photon hybrid system with microwaves and record the reflection spectrum 
S11. We describe the hybrid system in the Tavis-Cummings model with the 
Holstein-Primakoff transformation by two coupled harmonic oscillators, Eqs. (2.24), 
(2.26), and (2.31), 


Asys = hora a + hmm + hig (am? dE atm) ; (2.38) 
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that is used to derive an analytic model for the measured reflection spectrum of 
the system. 


While the Tavis-Cummings model does include the magnon and photon excitation 
number operators, we work in the first part of this thesis in a purely classical 
regime, where such a quantum description of the cavity magnon polariton would 
not be necessary. The model becomes important in Chaps. 4 and 8, where the 
superconducting qubit as the non-linear circuit element is introduced. 


A second, more intuitive approach to magnon-photon coupling treats the system as 
two coupled harmonic oscillators similar to classical physics lectures. Here, we find 
an eigenvalue equation of two coupled harmonic oscillators in matrix form [51] 


(“ Fa, soe : ta = (2.39) 
g w — Om X 


with constant coupling strength g and the frequencies of the photon and magnon, 
Õr and Wm respectively, which include complex damping terms. We make a classical 
rotating wave approximation, where we neglect higher order terms in w and find 
the eigenfrequencies of the coupled system in the vicinity of zero detuning 


oa () at sy (r= Gm MY y 9 bit) 


We fit the measured dressed frequencies to obtain the bare cavity frequency and the 
current dependent magnon frequencies as well as the coupling strength. 


There is also a third description model of the magnon-photon coupling based on 
the electromagnetic viewpoint, which reveals best the underlying physics and is 
the basis of phase coherent effects such as spin-pumping [99], dynamic control 
of the coupling strenght [100, 101] or coherent time control [102]. However, this 
treatment of the coupled system is not relevant for the experiments performed during 
this thesis. 


2.3.1 Derivation of the coupled reflection spectrum 


We connect the system by coaxial cables to the outside world in order to probe the it 
with microwaves and measure the reflection spectrum. The probing of the system 
Hamiltonian is modeled by a bath of photons in the microwave feed line and an 
interaction term of the bath with the cavity photons where we consider only one 
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single mode within the cavity and no direct coupling between the bath photons and 
the magnon excitations (see Fig. 2.4), 


H= Psys T Apath + At (2.41) 
i= l dwħwbt bw, (2.42) 
Ain, =-ih / deco (bta — atu). (2.43) 


m Cae i. 


Figure 2.4: Hybrid system probing scheme The coupled cavity photon - magnon system is probed in 
reflection. The reflection spectrum is modeled with a continuous bath of photons in the feed line that 
couple only to the cavity photon mode. The excitation is then swapped to the magnon via dipolar 
exchange coupling. 


The interaction with the external bath is purely imaginary, since it acts as an 
additional source of dissipation for the cavity photons and has a negative sign to 
ensure a Hermitian Hamiltonian [51]. We employ the input-output formalism [103] 
and find the field operators for the incoming and outgoing probe field to calculate 
the reflection spectrum 


A 


Si = Pout | (2.44) 
bin 


From the Heisenberg equations of motion for the field operators 


db. i 


Se = 5 Ba, Al = —iwbu + cwû (2.45) 
dâ ite ip. y 2 A 
T = -> lâ, A] = -> (a, Hays] = o dw Cabe (2.46) 
a i a 
SE = -> Im, A] = = [m, Boys] (2.47) 
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we find by solving Eqs. (2.45) and (2.45) for a time tọ in the past 


A A y t É / 
bur =bu (to) eiet) + | dt! ceil") (t) (2.48) 
0 
då ie o ee 
Tr == = lâ, Asys] — sl dw cwe Wett), (to) 
00 t A A 
2 / dw / dt! | ewe-Na(r). (2.49) 
—00 to 


The last term in Eq. (2.49) describes the coupling from the bath into the cavity that is 
given by the corresponding linewidth xe. Using Fermi’s golden rule we identify 


2K. (Wy) = 27 [aw [cwl ô (wr — w). (2.50) 


Employing the Markov approximation, we assume a constant coupling over a 
frequency range close to resonance and find cy 7 c = yKe/27. The remaining 
integrals give a ô function and a factor 2m /2. We identify the second term in Eq. (2.49) 
with the incoming photon mode [104-106] 


bin (t) = 2 E dw e™™® (t-t) b. (to). (2.51) 


Choosing different boundary conditions for the time integral in Eq. (2.48) and 
defining an output mode with opposite sign compared to Eq. (2.51) we find two 
solutions for the equation of motion of the cavity operator 


= = > [â, Asys] + V2Kcbin — Kea (2.52) 
= = si fâ, Bsys] — V2Kebout + Kâ, (2.53) 

and the relation of the incoming and outgoing field operators 
Bin (£) + bout (t) = V2xcA (t). (2.54) 


This gives a set of two equations for the photon operator, Eq. (2.52), and magnon 
operator Eq. (2.47), 


a = (-ior =) a(t) ight (t) + V2Xebin (t) (2.55) 
a ON (2.56) 


We Fourier transform the operators into frequency space, leaving a factor of —iw on 
both left sides and solve the set of equations for á (w), 


â (w) = bin (w). (2.57) 
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With the input-output relation in Eq. (2.54) we find the reflection spectrum of the 
coupled harmonic oscillator system 


511 (w) = + =-14 (2.58) 
bin (w) i (wr, — wW) + Ke + TE 


Analog to Sec. 2.1.2 we introduce internal losses in the cavity and magnetic system 
by a complex renormalization of the resonance frequencies, wr — wr — ik; and 
Wm —> Wm — iKm [81] and find the reflection spectrum in the general case 


2 
Be = —. (2.59) 


i (wr — w) +K + 


i(wm—w)+Km 


2.4 Yttrium iron garnet 


We choose yttrium iron garnet (Y3Fe5O 2, YIG) as the magnetic material for all of 
our experiments. YIG, an insulating ferrimagnet, is the material of choice in many 
magnetism focused studies for its unparalleled low damping qualities and high 
Curie temperature of Tc = 560K. 


b 


Figure 2.5: Crystal structure and photograph of YIG (a) Conventional unit cell of YIG, where black 
spheres are yttrium, red spheres are oxygen, and green and blue are iron atoms on two different lattice 
sites (tetrahedral and octahedral, respectively). From [107]. (b) Photograph of one of our used YIG 
sample. The white ring is an artifact, showing the reflection the lighting equipment. 


While magnon lifetimes in pure iron or the alloy permalloy (Nig1Fe19) are on the 
order of nanoseconds [108, 109] and a spin-wave mean free path typically below 
10 um [110, 111], the low damping in YIG result in a lifetime of a few hundred 
nanoseconds and spin-wave propagation over centimeter distance. 
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These properties have led to a broad usage in research studying microwave 
magnetic dynamics [112], general wave dynamics [113], room temperature Bose- 
Einstein condensation [44], magnon logic elements [114], as well as cavity-based 
QED [115] optomagnonics [57], magnomechanics [116], and temperature dependent 
magnetization dynamics [117]. The low damping properties also resulted in the 
development of many analog microwave components such as filters, oscillators, and 
circulators based on YIG. YIG has a bcc lattice structure with 80 atoms per unit cell 
with a lattice constant of 12.4 A. The structure is depicted in Fig. 2.5(a) and a zoomed 
photograph of our sample is shown in Fig. 2.5(b). The magnetic properties originate 
in the two different sites of the Fe** ions which, differin the configurations of the 
neighboring oxygen ions. This results in a strong magnetic exchange interaction, a 
magnetic moment of a unit cell of 10 ug at low temperatures [118], and an effective 
spin number s = 5/2 [53, 115]. Due to its complexity a lot of research has been spent 
on the crystal growth process leading to well established single crystal processes 
since YIG cannot be found naturally. 


2.5 Magnetic losses 


Employing magnons in hybrid quantum circuits requires their excitation lifetime to 
exceed the exchange time between system components. The excitation lifetime is 
represented by the magnon resonance linewidth, which quantifies excitation losses. 
Measuring excitation losses in magnetic materials has become a standard procedure 
for material characterization by comparing the magnetic resonance linewidth. 
Depending on the physics language, this is expressed either by measurements of ôH, 
the linewidth of the magnetic resonance field at constant measurement frequency 
or by measurements of dw, the linewidth of the resonance frequency at constant 
magnetic field. The data can be converted between the two representations using 


ôw 
m dwy/dHy í 
with the differential dwy /0Hy representing the magnetic dispersion relation [119]. 


(2.60) 


We choose the language of microwave engineering and use the resonance frequency 
linewidth dw as the magnetic linewidth in this thesis, to which we refer to as 
Km. After the discovery of YIG in the late 1950s there have been multiple studies 
regarding the main underlying physical contribution to magnetic linewidth in 
different temperature regimes to further optimize the material quality. The most 
prominent contribution are found in material impurities, with other rare-earth (RE) 
elements replacing yttrium on a sample site, and scattering on the rough sample 
surface. At millikelvin temperatures, quantum defects in the form of two-level 
systems are found in YIG. Detailed experimental data is given in Chap. 7. 
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2.5.1 Slow-relaxing impurities 


Below room temperature, the linewidthin YIG increases with decreasing temperature, 
showing a peak between 40 K and 150 K [120-122]. This increase is explained by 
scattering at RE impurities within the YIG sample, the so called slow-relaxing 
impurities [119, 123, 124]. The two Fe** sublattices are strongly coupled and can 
be considered as one effective sublattice. RE ions that substitute on the yttrium 
sites couple to this effective sublattice with an exchange interaction about an order 
of magnitude smaller than the Fe*+-Fe** exchange interaction [122]. The RE-RE 
coupling can be neglected for small impurity concentration. The effective Fet 
sublattice couples to the non-interacting rare-earth ions individually. This coupling 
acts as a loss channel to the lattice for the uniform magnetic mode at a frequency w, 
broadening the resonance linewidth. 


The motion of the magnetization around the applied static magnetic field changes 
the exchange field acting on the parasitic RE ions. The time-dependent component 
parallel to the RE ion spin modulates their Zeeman energy splitting from the 
unperturbed frequency Q, the thermal equilibrium population changes, and the 
system relaxes back to equilibrium by energy level transitions. Due to the finite 
impurity relaxation time there is a time delay between transition frequency and 
equilibrium population, which results in the magnetic excitation loss. We model 
the anisotropic exchange by assuming the impurities with only two energy levels, 
which approximates them as effective spin 1/2 systems, with all impurity splittings 
at the same transition frequency Q with relaxation rate t. The linewidth due to 
slow-relaxing impurities «rg is found to be 


1 exp (AQ/kgT) WT 
T [exp (1Q/kgT) +1]? P +?" 


KRE X (2.61) 


with a typical impurity frequency on the order of (2/21. ~ 2 THz [119]. Considering 
a thermal equilibrium population of the impurities 


1 
gi exp (ħw/kgT) +1’ 


(2.62) 


the effect of slow-relaxing RE impurities freezes out towards single digit 
kelvin temperatures. 


2.5.2 Two-magnon scattering 


The dominant relaxation process for the magnetic excitation of the uniform mode 
going towards liquid helium temperatures is found to be two-magnon scattering 
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processes in imperfections on the sample surface [125-127]. The "two-magnon" 
expression originates in the second quantization formalism, where the process can 
be modeled as a magnon in the uniform mode being destroyed and a k 4 0 mangon 
at the same frequency being created [128]. Because the number of magnons stays 
constant, so does the z-component of the macro spin, but the destruction of a uniform 
mode magnon reduces the total spin by one. This makes it a scattering process 
only for the transversal component of the magnetic moment. The pit-scattering 
theory of Sparks, Loudon, and Kittel considers scattering at spherical diamagnetic 
inclusions at the sample surface. The linewidth contribution due to surface scattering 
Ksurface for a spherical sample of radius r with pits of radius rp, and the saturation 
magnetization Ms is well approximated by 


n? Y pit 


Ksurface % 2 Ms (2.63) 


r 


for the uniform precession mode [129]. 


2.5.3 Two-level systems 


At temperatures reaching the millikelvin regime the dominant loss mechanism in 
amorphous solids are two-level systems (TLSs). This rather generic loss mechanism 
was first described in experiments at ultrasonic properties of glasses [130-132], 
but has also been successfully applied to describe decoherence and fluctuations 
in superconducting quantum circuits [133, 134]. Additionally, TLS effects were 
observed in micro- and nanomechanical resonators [135-137], within amorphous 
reflective coating in optical devices [138] or field effect transistors [139]. 


In the standard tunneling model a two-level system is described as a double-well 
potential with a small energy difference ô and a tunneling rate 69, sketched in Fig. 2.6. 
The origin of this loss mechanism is not based on impurities but rather on the 
amorphous structure of the material. The TLS exists in one of the states and can 
change its state either by thermal activation or quantum tunneling. The general TLS 


Hamiltonian 
a 1 —ô ĉo 
H = 5 do 6 ) (2.64) 


is written in the basis of its eigenstates with energy E = y $2 + 65 and the Pauli 


Fe E S 
matrix 0, = 0 1 as 


Ĥ = -Eô;. (2.65) 
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The excitation loss in a TLS can be understood as a scattering process. If excitation 
and TLS are close in energy there is a resonant absorption process that depends on 
the occupation of the TLS. The general temperature dependence of the resonance 
broadening due to TLSs can be obtained by looking at the occupation number 
difference of the two states. 


energy 


Figure 2.6: Double-well potential modeling a TLS The energy of the two states is given by the ground 
state energy difference ô and the tunneling rate ôo. The barrier can be overcome by thermal excitation 
or quantum tunneling. 


With the total density of TLSs nyrs, the individual densities of TLSs in the lower 
and upper state, nı and nz, and assuming a Boltzmann distributed occupation ratio 
between nı and state nz, we find 


NTLS = n1 +12 (2.66) 

= exp (—E/kpT). (2.67) 

The TLS broadened resonance linewidth is proportional to the occupation num- 
ber [131, 140], 

KTLS X ÔN = M2 — N1. (2.68) 


With Eqs. (2.66) and (2.67) we find 


E 


For excitation frequencies w in the single digit GHz range the TLSs saturate at 
temperatures above a few kelvin, hw = kgT. The TLS occupation is not only 
dependent on temperature but also on excitation power P. The Hamiltonian in 
Eq. (2.65) is expanded by an interacting Hamiltonian modeling a drive field with 
transversal and longitudinal components. This system is equivalent to a spin 1/2 
particle in a magnetic field 


Ê ==hy (Bo + Br) -S, with ES = 7S x (Bo dj Br) (2.70) 
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with a gyromagnetic ratio y. This gives the well known Bloch equations [141] for 
the expectation values of the individual spin components. Calculating the magnetic 
susceptibilities using 

S; (t) x hyB"(t), (2.71) 


approximating a small amplitude in the interaction magnetic field, assuming only 
transverse coupling, and integrating the TLS energy with respect to the asymmetry, 
tunnel barrier and the transverse components we find an additional power dependent 
linewidth contribution [82, 142] 


KTLS X ==, (2.72) 


with the critical power P. at which the TLSs start to saturate. We obtain the total TLS 
linewidth contribution 


(2.73) 


where all material constants are combined into the not further specified linewidth x4 g. 


The microscopic origin of TLSs in any system is still not yet fully understood and 
part of active research regarding its nature and position in amorphous oxides. The 
possible theories range from a literal movement of single tunneling atoms or a group 
of atoms [143], dangling OH bonds [144], electron tunneling [145] to spins and 
magnetic impurities [146-148]. 
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The previous chapter describes a purely linear and classical hybrid system. In order 
to explore quantum features, we need an additional element in our system. This 
chapter introduces the needed nonlinearity, a superconducting transmon quantum 
bit (qubit). First, we give an introduction to the general concept of quantum bits and 
describe the concept of superconductivity and the Josephson effect, which makes 
our qubit nonlinear in the first place upon which our quantum bit is founded. After 
that, we elaborate the working principles and features of our superconducting qubit 
type, the transmon. The coupling of the quantum bit to the outside world via a 
microwave cavity photon concludes this chapter. 


3.1 The quantum bit 


A qubit is the quantum analogue to the classical binary bit, an artificial quantum 
mechanical two level system. The qubit is represented by its state vector |), a 
superposition of ground and excited states |g) and |e), respectively, 


|p) =a|g)+Ble), (3.1) 


with normalized complex mixing coefficients « and ß. The basis of the qubit 
eigenstates is spanned by |g) and |e), which can be represented by eigenvectors of 
the Pauli matrix 0, 


Is) = a) le) = a and ô; = E a . (3.2) 


With the qubit states quantized along the z-axis we find the creation and annihilation 


operators 04 = (ôx + id) /2, where ô- |e) = |g). 


Typically, the Bloch sphere representation is used to illustrate arbitrary qubit states, 
but also the impact of decoherence. Any pure qubit state |p) is represented by a 
point on the sphere surface characterized by a Bloch vector and corresponding Euler 
angles 0 and y, 


|y) = cos : |g) + sin seit le), (3.3) 
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depicted in Fig. 3.1. In analogy to Sec. 2.5.3 we attribute a transition energy hw, 
to the qubit states. The time evolution operator rotates a pure qubit state with the 
qubit frequency wq around the z-axis, resulting in a Larmor precession. Usually this 
rotation is accounted for by looking at the qubit Bloch vector from a rotating frame. 


Figure 3.1: Qubit state on a Bloch sphere A pure qubit state |p) is represented by a point on the Bloch 
sphere's surface. The ground and excited state of the qubit are eigenvectors of 0. By convention, the 
ground state is located on the Bloch sphere’s positive pole. The Bloch vector of the qubit state is defined 
by the Euler angles 6 and 9. 


The qubit dynamics in the presence of energy relaxation and decoherene are described 
by the Lindblad master equation [96, 149] that results in the time evolution of a 
qubit state driven by a frequency wp. The two relevant decoherence mechanisms 
are the energy relaxation and pure dephasing. Energy relaxation acts on the qubit 
with the decay operator 1; = 6_, relaxing any excited state to the ground state at 
a rate I]. Pure dephasing results in phase fluctuations of the Lamor precession in 
the equatiorial plane at a rate Ir. Here, the decay operator is given by Î = 6,/ V2. 
In the Bloch sphere representation, the energy relaxation operator ö_ maps the 
|e)-component of the state to the ground state eigenvector, while the dephasing 
operator describes a Larmor-type precession of the state around the quantization 
axis, describing the dynamics in the equatorial plane. 


Solving the Lindblad equation [80] we find two characteristic time scales for the 
decoherences. Tš = 1/I¥ describes the timescale on which the qubit loses all 
phase information and T; = 1/1; is the energy relaxation time.. Neglecting pure 
dephasing, I; = 0, we find the dephasing time limited by 


MEM (3.4) 
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The dephasing is measured by performing a Ramsey experiment [150], where 
the qubit is prepared in a state on the equatorial plane and rotated back on the 
quantization axis before reading out its state. Preparing the qubit state along the 
x-axis, the probability to find the qubit in its excited state after the second rotation 
decays exponentially with the dephasing rate If = IT 1 + I+, modulated by cos Awt 
with Aw = wq — wp with the drive frequency wp. The oscillation due to the frequency 
detuning Aw # 0 are referred to as Ramsey fringes. 


Driving the qubit by an oscillating drive or pulse, its state vector rotates continuously 
between |g) and |e). The qubit drive is modeled by 


Ha = hOR cos (wpt) Ox, (3.5) 


with the Rabi frequency Qpr [151]. Rotating the qubit state continuously allows us 
to calibrate pulse lengths required to prepare a specific state. Transforming the 
Hamiltonian into a frame rotating with the drive frequency and applying the rotating 
wave approximation by neglecting fast rotating terms, we find the probability to 
find the qubit in the ground state using 


Pie = (sl |a) (3.6) 


with the time evolution operator in the rotating frame Ü (t). For zero detuning we 
find the ground state probability [152] 


Pig) (Aw = 0) = > (1 + cos Ort). (3.7) 


The Rabi frequency increases with increasing drive power [153]. 


3.2 Superconductivity and Josephson junction 


The transition from the concept and theoretical description of a qubit towards 
fabricating a measurable and controllable device that exhibits quantum coherence 
requires a manageable amount of effective degrees of freedom. We chose a solid 
state approach by using superconductors as our engineering platform for realizing 
a qubit. Superconductors essentially allow for a macroscopic quantum coherent 
flow of charge due to a long-range order parameter [154]. Although the typical 
implementation of a superconducting qubit has a spatial size on the mm length 
scale containing around 101? atoms it is still possible to describe them with a single 
degree of freedom [79]. For the scope of this thesis it is sufficient to focus on the 
macroscopic phenomena of superconductors [155] without exploring the details 
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of the underlying microscopic description by the BCS theory, whose developer 
ultimately were awarded a Nobel prize in 1972 [156]. 


At a critical temperature T. a superconductor undergoes a second-order phase 
transition where its electrons forms a Cooper pair condensate. Cooper pairs consist 
of two electrons with opposite spin and momentum, resulting in zero effective spin 
and a boson-like nature. The weak net attraction between the two usually repulsive 
electrons is mediated by an interaction with lattice phonons. The Cooper pairs 
undergo a Bose-Einstein condensation, leaving them in a macroscopic quantum 
state described by a single macroscopic wave function Y (r) = |'¥| e!?) [85] with 
the collective phase ọ (r) and the amplitude || being related to the Cooper pair 
density. The flowing supercurrent is determined by the spatial dependence of y 
and is a quantum coherent motion of the condensate as a whole. The Cooper pair 
condensate is broken by a finite pair breaking excitation 2A per Cooper pair as lowest 
possible excitation. Scattering processes below this minimal energy are prohibited, 
resulting in a non-dissipative current flox. This behavior was first observed in 1911 
in mercury [157]. The scale of the gap is on the order of tens of GHz or several kelvin, 
i.e. for aluminium, Aa] = 41 GHz [85]. In addition to a vanishing electrical resistance 
a superconductor exhibits perfect diamagnetism [158], expelling magnetic fields 
from the bulk. 


A strip of superconductor forms a harmonic LC oscillator similar to a metallic 
antenna. The resonance frequency (and higher harmonics) is determined by the 
spatial extension and the effective dielectric constant of the medium. The inductance 
and capacitance of such a system are given by a geometric part, storing energy 
in magnetic and electric fields and a kinetic part, storing energy in the motional 
degrees of freedom of the charge carriers. Frequent electron scattering in normal 
metals prevents the storage of any significant amount of kinetic energy, but due to 
the suppressed low-energy states in superconductors, the kinetic inductance term 
becomes more important, especially in non-uniform superconducting materials 
as oxidized aluminium [159-161]. Similar to the three dimensional cavity case in 
Sec. 2.1 we can couple superconducting resonators via their electromagnetic dipole 
moment to the outside world. 


Connecting two superconductors by an insulating tunnel junction (see in Fig. 3.2(a)) 
leads to an overlap and hybridization of the two individual macroscopic wave 
functions due to tunneling effects. This ideal tunnel junction, together with the 
intrinsic capacitances of the electrodes, creates a Josephson junction [162] as depicted 
in Fig. 3.2(c). The Schrödinger equation for the hybridized wave function connects 
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the supercurrent across the barrier I; with the superconducting phase difference 
$ = $1 — Q2 and the critical supercurrent Ke, 


I = ksing (3.8) 


and connects the voltage drop across the junction V with the flux quantum Dy = h/2e, 


dd 
V=> 319 (3.9) 
a b c 
Epot 
i 
|W] exp (ime (r 2Ej |1) 


Figure 3.2: Josephson junction and anharmonic potential (a) Two superconducting electrodes connected 
by an insulating tunnel barrier form a Josephson junction. The current and voltage dynamics of the 
tunnel current are determined by the phase difference ¢ = 91 — 2 leading to a nonlinear inductance. 
The amplitudes of both macroscopic wave functions are considered identical. (b) The non-linear 
inductance leads to an anharmonic potential and non-equidistant eigenenergies of the corresponding 
states of the oscillator. The energetically lowest states are used as the two levels of a superconducting 
qubit. A sketched harmonic potential shows the impact of the nonlinearity on the potential landscape. 
(c) Circuit diagram of a Josephson junction. The cross symbolizes the insulating tunnel junction, with 
the intrinsic self capacitance in parallel. 


The Josephson junction allows for a zero-voltage supercurrent of Copper pairs to 
tunnel across the barrier, before the junction becomes resistive at a pair-breaking 
current Ie. Inserting the Josephson equations Eqs. (3.8) and (3.9) into the definition 


of the response of a inductor V = —L4 we obtain a non-linear inductance for the 
Josephson junction, 
Do 
L = — 1 
1 (9) 2T Ie cos p (ao 


The corresponding Josephson energy is given by 
E = >k (3.11) 


and we find a charging energy Ec = e?/2C, accounting for the isolated superconduc- 
tors acting as a small capacitor. The sinusoidal inductance alters the total inductance 
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of the equivalent LC circuit and the energy levels of the corresponding oscillator 
become non-equidistant as depicted in Fig. 3.2(b). This anharmonic oscillator is the 
basis for superconducting qubits, where the lowest two energy levels are used as the 
ground and excited states of a qubit. 


3.3 The transmon qubit 


In this thesis we use a transmon qubit as the qubit type of choice. The transmon 
qubit [163] builds on the design of the single Cooper pair box (CPB) [14, 164] which 
consists of a superconducting island that is connected to a Cooper pair reservoir via 
a Josephson junction. Both, CPB and transmon qubit are described by the excess 
Cooper pairs N on the island and the superconducting phase ¢ across the Josephson 
junction. The number operator is given by 


Ñ=} NIN) (N| (3.12) 
N 


and its conjugate phase operator ¢, 


[N, $] = -i. (3.13) 
The basis states are related via Fourier transformation 
1 pom , 
N) = — ING 
IN) = 5 [ape iN? |p) (3.14) 
le) = Y el? IN), (3.15) 
N=-—00 


Cooper pair tunneling through the Josephson junction is described by [165] 

A Er E 

Ay=-2 Y IN (N+1|+|N +1) (NI, (3.16) 
N=-o 


where N is the number of Cooper pairs on one electrode connected to the Junction, 
later the island electrode. This allows to define a charge operator 


Ó = —2eN. (3.17) 


This results in the Josephson junction Hamiltonian from Eq. (3.16) to become 


f Er 1 27 A : 
A=-d f ao (er +e) 0) 0) 
= -Ejcos® (3.18) 
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For the intrinsic Josephson junction capacitance we find 


Ac = — =4EcN’. (3.19) 
This results in an effective Hamiltonian for the Josephson junction 
Ay = Ay T Ae = 4EcN? — Ej cos $. (3.20) 


Transitioning to the CPB we introduce an additional term representing the offset 
charge Ng on the superconducting island and find a Hamiltonian similar to the single 


Josephson junction, 
A A 2 A 
H=4Ec (N — ng) — E cos p. (3.21) 
a b c 
z E/Ec=0.5 E E/Ec=5 E/Ec=50 
2 
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Figure 3.3: Qubit eigenenergy dispersions for different ratios E; /Ec Eigenenergies E,, for m = 0, 1, 2 
as a function of the offset charge ng for different characteristic energy ratios Ej /Ec. The energies are 
normalized to the first transition energy Eo; at the degeneracy point ng = 1/2. (a) In the charge regime 
of the Cooper pair box, the transistion frequency is sensitive to changes in the offset charge, with 
charge sweet spots at half-integer ng. In the transmon case (c), the charge dispersion flattens and the 
transition frequencies become insensitive to the offset charge. This comes with the cost of a reduced 
anharmonicity, described by a weak anharmonic potential. (b) shows the transition regime. Figure 
similar to [163]. 


Figure 3.3(a) shows color coded the lowest three eigenenergies E; of the CPB. We 
treat the lowest two states as our qubit, with the qubit transition frequency being the 
energy difference between them. Usually, the CPB is biased to operate at half integer 
Ng, creating a coherent superposition of charge states. 


With the parabolic energy dependence on the offset charge, the qubit frequency is 
sensitive to charge fluctuations. The transmon qubit was proposed to operate in a 
regime where the ratio of Josephson energy and charging energy is increased [163], 
E;/Ec 2 50. This increased ratio reduces the anharmonicity (non-equidistant energy 
level spacing) but decreases the charge noise sensitivity. Figure 3.3 shows the 
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eigenenergies for the CPB Hamiltonian for different E; / Ec ratios. For large ratios the 
charge dispersion is nearly flat and the qubit transition frequencies are effectively 
constant with respect to the offset charge. A transmon qubit design eliminates the 
need for a gate to tune the qubit into its charge sweet spot at half integer ng to 
minimize noise effects. The large characteristic energy ratio is realized by adding 
a large shunt capacitance in parallel to the Josephson junction that increases the 
charging energy with Ec «x C -1 The CPB Hamiltonian, Eq. (3.21), is still valid, but 
with the dominant Josephson term the phase instead of the charge becomes the 
well-defined variable due to the uncertainty principle. 


With the nearly harmonic energy level structure, the transmon resembles a simple 
harmonic LC-resonator where the inductance is replaced by the nonlinear Josephson 
junction. This nonlinearity creates a weakly anharmonic potential and allows the 
distinct excitation of the individual transitions. The transmon eigenstates and 
eigenenergies are approximated by expanding the cosine term in the potential 
energy up to fourth order in $, which is valid for large E; /Ec. The system resem- 
bles a harmonic oscillator with a quartic perturbation that describes the leading 
order anharmonicity, 


A = Ech (6% | z) EE (64697, (3.22) 


A 


with the approximated harmonic oscillator creation and annihilation operators, b 


and bt. The eigenenergies are corrected in the leading quartic order by 


A iy aries O BO Oy E 
EN = — SE (l(b +54) = —=E (6? +6j +3), (3.23) 


where |j) are the purely harmonic oscillator states. The transmon eigenenergies of 
state m are approximated by 


1\ E 
Em ~ —Ey + ,/8EcE) (m | 5) 5 (6m? 1 6m 4 3), (3.24) 


with the transition frequencies Emn = Em — En. As with the CPB, the lowest two 


states are used as our qubit two-level system. Overall, the transmon can be viewed as 
an artificial atom with multiple individual level transitions. The absolute transmon 
anharmonicity is given by the difference between the first two energy levels and is 
approximated by the charging energy, 


a= E12 = Eoı x= —Ec, (3.25) 


with the relative anharmonicity 


& = —4/8E¡/ Ec. (3.26) 
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3.4 Jaynes-Cummings model 


For manipulation and readout, we couple the qubit via photons in the form of 
confined microwave cavity modes to the outside world. In general, this is described as 
light-matter interaction of an atom and a single resonator mode. In our experiments 
this is realized via the electric dipole of the transmon qubit, à « ôy, which is in 
parallel to the electric field of the cavity, Ê œ ât + â. We find a coupling Hamiltonian 


A. = hg0, (a + a) (3.27) 


that describes the transverse coupling between atom and resonator, with ôy having 
only off-diagonal terms. Approximating the atom to only two levels with energy 
difference hwg that interact with a single quantized mode of a microwave cavity at 
frequency wr we find the quantum Rabi Hamiltonian [151, 166] 


> 1 
A/h = sword: + wrt" + gör (at +a), (3.28) 


with creators/annihilators in the Fock space of the oscillator mode ât/â. The 
systems is usually described in a regime, where the coupling is small compared 
to the frequencies, g < wr, Wo. We use ô+ = 1/2 (0, + 0-_) and apply a rotating 
wave approximation to Eq. (3.28) where we neglect the counter rotating terms 


proportional to ô+ ât and 0_4. We remain with the Jaynes-Cummings Hamiltonian 
for the interaction of a qubit with a microwave cavity mode [167], 


; 1 
Pjc/ħ = swe: + writ 4 g (ate +644). (3.29) 


3.5 Dispersive qubit state readout 


The transverse coupling of qubit and readout cavity enables us to determine the 
qubit state by meseuring the cavity frequency. Having a large detuning between 
qubit and readout resonator frequency A = wo — Wr compared to the coupling 
strength, g/A < 1, we are able to perform a dispersive readout of the qubit state 
[16, 17]. In the dispersive regime we find a shift of qubit and resonator frequency by 
a factor x = g?/A for a perfect two level system, which allows for the qubit state to 
be measured along its quantization axis in a quantum nondemolition measurement 


[168]. Due to the low anharmonicity of the transmon, its dispersive shift consists 
2 E c 


of two contribution with opposite signs, x = Xo1 — X12/2 = TATA for large 
E; / Ec [163]. Applying a unitary transformation 

aea are 

U= exp (E (aâ 9-) } (3.30) 
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diagonalizes Eq. (3.29). Expanding it up to the second order of x we find [17] 
yt Wor» l ve) ata 
UAycU" /h = aa (w. + x6z) aA, (3.31) 


with the renormalized frequencies w}, = woı + X01 and wi = wr — X12/2 due to the 
cavity qubit interaction. The dressed frequency of the readout resonator depends 
on the state of the qubit, w(0,) = wr + x. By probing the resonator, the qubit 
Bloch vector collapses into one of the eigenstates of the ô; operator, |0) or |1). This 
interaction with the quantum mechanical system minimizes this impact on its time 
evolution and preserves the eigenstates of a qubit. Such a measurement scheme is 
referred to as a quantum nondemolition measurement. 


We plot the transmon population on a scale from 0 to 1, with 0 corresponding to the 
the qubit being in the ground state |0) and accordingly for 1, and repeat the single 
projective measurements for statistic purpose. 
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In this chapter we describe the combination of the physical systems discussed in 
Ch. 2, coupling of magnons and photons, and Ch. 3, coupling of a transmon qubit, 
an artificial atom, to photons. With the qubit as a non-linear element in the system 
and using a photon mode in a microwave cavity as a bus system, we are able to 
create a nonlinear hybrid system containing magnetic excitations. 


We explore the quantum hybrid system in two configurations. First the dispersive 
regime, where qubit and magnon are on resonance, far detuned from the cavity [65, 
66]. Here, the effective qubit-magnon interaction is mediated by an effective coupling 
via the resonator coupling mode. Such a system with an effective qubit-harmonic 
oscillator coupling was studied in the past using photons regarding photon number 
states [169, 170], coherent superposition states [171] and Schródinger cat states 
[172]. Recently, the first experiment regarding counting excitation number states 
was successfully carried out using magnons [67]. In the second configuration we 
describe a strongly coupled tripartite system, where all modes are degenerate and 
the strong coupling is directly observable in the reflection spectrum. 


We discuss the overall coupling scheme in the dispersive regime regarding the 
expected effective coupling strength between qubit and magnon ¿qm and describe 
the eigenfrequencies of the tripartite system by a set of three coupled harmonic 
oscillators and its reflection spectrum by expanding the input-output formalism to 
three participating components, where the transmon energy spectrum is truncated 
to only include the ground and excited qubit state. The chapter concludes with a 
short discussion of the influence of the applied magnetic field to the qubit frequency. 


4.1 Dispersive regime 


In the dispersive regime, the qubit frequency and the Kittel mode frequency are far 
detuned from the cavity resonator. Since there is negligible direct qubit-magnon 
coupling, the effective coupling is mediated by a virtual photon coupling using 
the cavity modes, with the nearest mode bearing the greatest share. The system is 
probed and read out similar to the bare qubit in the dispersive regime, monitoring 
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the state dependent resonance frequency of a higher cavity mode. On qubit-magnon 
resonance and given a strong effective coupling, the dispersive shift shows a magnon- 
vacuum-induced Rabi splitting [65]. For the system Hamiltonian we model the TE10p 
mode resonance of the cavity and the Kittel mode resonance as harmonic oscillators 
and treat the qubit as an anharmonic oscillator with anharmonicity « ~ —Ec < 0, 


A 2 
Hays /h = Y wropapay + (0 5) 49+5 (9) | + omir, (4.1) 
p 


with the resonance frequency @10p of the TEjop cavity mode and the qubit and 
Kittel mode frequencies wg and wm, and the corresponding creation/annihilation 
operators ât /ap, G ĝt /g, and mt /í1 for the cavity mode, the qubit, and the Kittel 
mode En. respectively. We neglect direct qubit-magnon coupling and write the 
interaction Hamiltonian for the remaining components, 


Anı/h = 280107 Gi dp + ata) + 28m, 10p (m! Ap + aj) , (4.2) 


with gq10p and Sm,10p being the coupling strength of the TEj9p mode with the qubit 
and the Kittel mode, respectively. 


In the dispersive regime, the detuning of qubit and Kittel mode from the cavity 
modes much larger than any other energy regarding coupling strengths. Tun- 
ing the Kittel mode frequency towards the qubit transition, in the regime where 
8q10p, Sq lop < |W — Wm|, the interaction Hamiltonian is expressed in the rotat- 
ing frame of the system Hamiltonian under perturbative approximation to the 
first order [173] as 


p I>0 


apm ñ, (4.3) 


where the first two terms describe dispersive shifts due to the cavity-magnon 
coupling, the third and fourth term describes the qubit frequency’s Lamb shift and 
ac Stark shift, and the last term describes the static interaction between the cavity 
modes and the Kittel mode [66]. 


When the qubit and the Kittel mode frequency are degenerate, the detuning is 
smaller than the couling strengths and their coupling is mediated through virtual 
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photons in the far detuned cavity modes. The qubit-magnon detuning smaller 
than the coupling strength allows to transform the interaction Hamiltonian in the 
corresponding rotating frame using the Schrieffer-Wolff transformation [47], which 
eliminates the virtual photons. The effective qubit-magnon coupling on resonance is 
approximated by, 


Pim /t = Sqm (M0 + Om) , (4.4) 
with an effective coupling strength 


e 8m,10p8q,10 
Sqm = | Ta, (4.5) 
p W10p = wq 


4.2 Tripartite system: light-matter-spin 


In contrast to the dispersive regime, the detunings of the component frequencies in 
the tripartite system are smaller than the coupling strengths. The system is inspired by 
strong light-matter interaction experiments [16], where coherent interaction between 
qubit and microwave photon create an anticrossing in the reflection spectrum due to 
vacuum Rabi oscillations. Additionally, we tune the magnon frequency through the 
dressed qubit-photon states to further split up the reflection spectrum. The general 
coupling scheme is depicted in Fig. 4.1. Increasing the coupling from the uncoupled 
case (subfigure (a)) to the full qubit-resonator coupling (subfigure (b)) creates the 
dressed states. The magnon adds a third degree of freedom into the system, resulting 
in three modes with field dependent detunings A; and A) (subfigure (c)). 


a b c 
uncoupled cavity- cavity-qubit cavity-qubit-magnon 

qubit system coupling coupling 

> 

> 

$ Ao (B) 
| | A; (B) 

Ig =0, B=0 Jg #0, B=0 Io #0, BAO 


Figure 4.1: Schematic energy diagram of the strongly coupled tripartite system With nearly degenerate 
photon (yellow) and qubit (blue) frequencies the system forms dressed mixed states in the strong 
coupling regime that are not tunable (b). (c) The magnon (red) frequency is tuned with the applied 
magnetic field and the coupled system has three degrees of freedom that form three hybrid states with 
variable frequency spacing A; (B) and A; (B). (a) shows the zero coupling case in the presence of a 
strong probe drive. 
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We expand the general coupling scheme between two harmonic oscillators, Eq. (2.39), 
to three components, with the qubit treated as a harmonic oscillator in the weak- 
excitation approximation, and find the 3 x 3 coupling matrix 


W_ Wy Srm Sqm 
A= Sm WOW gar p (4.6) 
Sqm gar w Wg 


with the eigenvalues resulting in the dressed frequencies w /2/3. With the 
characteristic polynomial being a cubic function we obtain the eigenvalues 
by numerical calculation. 


4.2.1 Reflection spectrum 


We describe the reflection spectrum S11 of the tripartite system using the input- 
output formalism as described in Sec. 2.3.1. For this, the system Hamiltonian in 
the rotating wave approximation is expanded by the qubit contribution, that is 
approximated with the 6, operator as a perfect two level system, and the respective 
coupling strengths, 


A hw 
Asys = hw ât + hwm m + == 


+ga (010 +0,4) 4 gem (mia + atin) +gqm (ite +65), (4.7) 


with the frequency and creator /annihilator operators of the resonator, qubit, and 
magnon, wr and a/a, wg and 64/6_, and wm and ió Im, respectively. The indexes 
in the different coupling strengths indicate the participating components. The 
Hamiltonians describing the photon bath in the feed line and the system-to-bath 
coupling are the same as in Eqs. (2.42) and (2.43), respectively, which results in the 
same equation of motion for the resonator operator as in Eq. (2.52). 


We use the weak-excitation approximation and assume the qubit to be in the ground 
state, (07) = —1 [81, 174, 175], and find the set of three differential equations 


E = ( iw, =) a(t) igrmift (t) — iggpd— (t) + V2kcbin (t) (4,8) 


dt 2 

dí : x 3 a a A 

T = —iwm/t (t) —igrmá (t) — iggmó- (t) (4.9) 
dö_ RE RER A z 
E = —iwgÎ- (t) — iggrá (t) — İgqm (t) E (4.10) 


The system is then Fourier-transformed and solved for the individual annihilation 
operators. We insert the resonator annihilator into the relation of the reflection 
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spectrum to the incoming and outgoing field operators bin /bout, Eq. (2.44), and 
include the complex renormalization of the resonance frequencies wj — wj — ix; to 
account for internal losses in component j, 


2Ke 


$1, = -1+ - , with 4.11 
! R | Sar i gan 2igqr8rm&qm ( ) 

Q+ sam eco MQ+ 85m 
R =1i(0r— 0) +n, (4.12) 
Q =i (w¿— w) + Kg, and (4.13) 
M=i(wm — w) + Km. (4.14) 


Comparing this to the spectrum containing two components, Eq. (2.59), we identify 
the second term in the denominator as resonator-qubit coupling with a subsequent 
coupling to the magnon. Similarly, the third term describes the resonator-magnon 
coupling with subsequent coupling to the qubit. The last term in the denominator is 
attributed to a correction of the qubit-magnon coupling that is counted in twice. 


The approximation to truncate the qubit to its two lowest level breaks down when 
other qubit levels are excited by multi-photon transitions or due to thermal excitation. 
In this case, the Hamiltonian in Eq. (4.7) with one qubit frequency is not valid anymore 
and the full transmon Hamiltonian, described in Sec. 3.3, with its anharmonic level 
structure and coupling strengths of the individual transitions to the other system 
components is used and we simulate the spectrum numerically. 


4.3 Superconducting qubit in magnetic field 


By tuning the frequency of the magnon excitation we expose the superconducting 
transmon qubit to an in-plane magnetic field. We work with a 3D transmon qubit 
design, see Sec. 5.2, with one single Josephson junction whose critical current I. 
is periodically suppressed by the magnetic field which results in the shape of a 
Fraunhofer diffraction pattern [176], 
( B 
sinc Be; ) 


with Bg, being the magnetic field corresponding to one flux quantum in the 
effective junction cross section area. Approximating the sinc function, sinc (x) == 


K (B) at : (4.15) 


1- = + O(x*), and with Ej « Ic, we assume a linear dependence with negative slope 
of the fundamental qubit transition frequency on the applied magnetic field. The 
single qubit fabrication step eliminates unwanted additional Josephson junctions, 
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whose large overlapp surface results in a higher sensitivity to applied magnetic 
fields and the I, dependence is described as a complex trigonometric relation of the 
two junctions [177]. The coherence of the transmon qubit decreases with increasing 
magnetic field, with the T; time being reduced by dissipation due to entering flux 
vortices in the aluminium thin film and their movement due to the rf current [178] 
and dissipation through quasiparticle losses [179-181]. 
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This chapter presents details about the sample fabrication process and the used 
experimental setup during this work. Our samples consists of three components, a 
commercially obtained magnetic sample, a machine-milled copper cavity resonator, 
and a nanostructured superconducting qubit. We briefly discuss the details of the 
magnetic sample, provide the basic fabrication concept for superconducting qubits 
using a Josephson tunnel junction, and discuss the cavity resonator. In order to access 
the quantum regime with our samples we perform the experiments at millikelvin 
temperatures. This strongly suppresses classical thermal noise and higher energy 
levels do not get populated thermally. We describe the cryogenic measurement setup 
and briefly discuss the working scheme of a *He/*He dilution refrigerator. The 
samples are probed by microwave excitations with respect to different magnetic bias 
fields, which requires certain hardware and wiring components. After that, we focus 
on the employed measurement technique using spectroscopic and time-resolved 
measurement schemes as well as the used measurement software package. 


5.1 Magneticsample 


Throughout this thesis we use two different commercially available magnetic YIG 
samples that we obtained from Ferrisphere, Inc.!. The single crystal YIG samples come 
diamond polished to a spherical shape with diameters of 0.5 mm and 1.0 mm, giving 
them an estimated 50 nm surface roughness. The spheres are already pre-mounted 
to a ceramic beryllium oxide rod oriented along the [110] crystal direction. 


5.2 Qubitsample 


The qubit samples fabricated during this thesis are in-house fabricated in the 
cleanroom facilities of the Nanostructure Service Laboratory of the Center of 


1 Ferrisphere, Inc., https: //www.ferrisphere.com - accessed Oct. 18, 2019 
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Functional Nanostructures at the KI'T. A detailed step-by-step list of the different 
cleanroom fabrication steps can be found in Appendix A. 


All samples are fabricated on intrinsic silicon wafers with a thickness of 380 um and 
a dielectric constant of e, = 11.5. Prior to any fabrication the substrate is cleaned 
with N-Ethyl-2-pyrrolidone (NEP) to remove organic residues. The Josephson tunnel 
junction and the remaining qubit structures are fabricated using electron beam 
lithography and shadow angle evaporation. 


The basic principle of the Dolan bridge technique [182] is to create a free-hanging 
bridge of resist and to evaporate metal from two different directions onto the substrate. 
This creates two electrodes with an overlap area forming the Josephson junction. We 
spin coat a 20 x 20 mm? wafer chip with a dual resist stack of LOR/PMMA. The 
resists have different critical doses regarding the electron exposure. At a certain 
dose, only one of the exposed resist is removed by the developer afterwards, which 
results in washing out the LOR under a free hanging PMMA bridge. 


The deposition of both metal layers is done thermally by evaporating an aluminium 
(Al) target in a PLASSYS Bestec MBE550S machine. Leftover resist residue is cleaned 
by a cleaning plasma with a 8 : 1 argon/oxygen mixture plasma before evaporating 
the metal. The the first Al layer is evaporated with a positive tilt angle before an 
insulating tunnel barrier is formed by dynamic in-situ oxidation. A controlled 
amount of oxygen is led into the chamber and the background pressure is adjusted 
via continuous oxygen input and pumping output. The 30 min exposure to the 
oxygen creates a ~ 1 nm thin insulating AlOx layer. The second Al layer is evaporated 
under a negative tilt angle and forms the Josephson junction, the connection between 
both electrodes. The two evaporation steps form the rest of the qubit as well. 


We evaporate 30nm and 50nm Al, resulting in an overall thickness of 80nm 
throughout the chip. The remaining resist is removed in a lift-off process, leaving 
only the written metallic structure on the chip. The chips are diced into 14 x 2 mm? 
strips, each containing a single qubit. 


5.3 Cavity resonator 


The cavity resonators are machine milled by our in-house mechanical workshop out 
of oxygen-free copper. For the experiments we use two similar cavities only differing 
in one dimension of the resonance volume. The nominal resonance volumes are 
83 x 40 x 3mm? and 83 x 37 x 3mmY, respectively, with rounded 1.5mm corner 
radii in the xz-plane due to fabrication restraints. The boxes feature a cut in 
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parallel to the xz-plane to access the cavity cutout and insert the other components 
of the sample. 


A hole with a diameter of 1 mm is drilled into one end of the box at half height to 
hold the magnetic sample in place. This is a point of maximum magnetic field of the 
resonance photon and the point where the external bias field is applied. The lower 
part of the opened box features a cutout to hold the qubit chip which is oriented 
in parallel to the electric field component of the resonance field. In order to lower 
the impact of the applied magnetic field to the performance of the superconducting 
qubit by worsen its coherence properties, the qubit chip and the magnetic sample 
are mounted on opposite ends of the resonator cutout. 


The cavity resonance gets excited by an 502 matched straight terminal SMA 
connector, with the inner conductor extending into the cavity volume in parallel to 
the electric field. The bottom part of the boxes feature several M3 threads to hold it 
in place in the cryostat and establish a thermal connection. The assembled sample is 
depicted in Fig. 5.1(c). 


5.4 Cryogenic setup 


All experiments are carried out in an Oxford Instruments Triton 200 dry dilution 
refrigerator, depicted in Fig. 5.1(a). 


The cryostat setup consists of five gold-coated copper plates at different temperatures. 
The two highest plates, the PT1 and PT2 stage at about 50 K and 3K, respectively, 
are cooled mechanically via a pulse tube cryocooler. Pressurized helium is used as a 
heat exchange gas and pumped from a room temperature compressor to the heat 
exchangers at these two plates. Releasing the pressure of the helium, the gas picks 
up heat at the heat exchangers and cools the connected PT1 and PT2 plate. The heat 
of the compressed gas is drawn off at room temperature. The three lower plates, 
still plate, intermediate plate, and the base plate containing the mixing chamber, are 
cooled via a secondary pre-cooler circulating the *He/*He mixture. The mixture is 
cooled at both pulse tube stages and picks up heat at the other plates on its way 
back to room temperature. We pump the volume inside the cryostat to vacuum to 
avoid any thermal contact between the different plates via air molecules. 


Cooling the fridge further down towards millikelvin temperatures, we employ a 
standard technique by using the phase separation of the He mixture. The warm 
mixture (75 L overall volume with nearly 15 % of 3He) is injected into the circulation 
lines with a pressure of about 2.5 bar. The mixture condenses in the still volume 
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making use of the Joule-Thomson effect, cooling the still, intermediate, and base plate 
down to about 1.2 K. The circulation process starts when all mixture is liquefied. 


Figure 5.1: Opened dilution refrigerator and sample box (a) Picture of the opened dilution refrigerator 
with the sample located at the very bottom. The temperatures of the stages decreases from top to 
bottom. (b) We apply an external magnetic field at the position of the YIG. A pure iron yoke with a 
superconducting coil and additional SmCo permanent magnets provide the magnetic field strength to 
tune the Kittel mode resonance frequency to the gigahertz regime. (c) The used copper box features a 
milled cavity, the YIG sphere mounted to a BeO rod positioned in an anti-node of the magnetic field, 
and the superconducting transmon qubit positioned in parallel near the maximum anti-node of the 
electric field of the cavity. The distance between YIG and qubit decreases the magnetic field influence 
on the superconducting Al. 
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Below a critical temperature of about 870 mK, the helium mixture undergoes a 
spontaneous phase transition in a 3He-rich and a *He-rich phase [183], which is 
forced out of equilibrium and the system extract energy from the environment to 
regain equilibrium. The *He-rich phase floats on top and consists of almost pure 
He, with the *He-rich phase containing about 6.6 % of “He. The circulated He is 
pumped into the *He-rich phase and disturbs the equilibrium concentration. The 
system relaxes back into equilibrium by bringing 3He into the ®He-rich phase, a 
process that effectively absorbs thermal energy from the base plate. The almost pure 
3He is pumped out of the still volume, cleaned at room temperature using a liquid 
nitrogen cold trap, and injected back into the mixing chamber from above. The 
warm gas is cooled on its way down by exchange cooler, using colder “He coming 
from the mixing chamber to the still. A dilution refrigerator reaches single digit mK 
temperatures employing this cooling process. Without additional sample wiring we 
are able to measure a base temperature of 9.8 mK in our cryostat using a %Co-based 
nuclear orientation primary thermometer, thereby calibrating our ruthenium oxide 
resistance thermometer for further cool-downs. Fully equipped our cryostat reaches 
working temperatures of about 25 mK. 


The cooper box with the cavity resonator and all other components of our sample 
are thermally and mechanically anchored with additional copper connections to the 
base plate of the cryostat. 


5.4.1 Magnetic field bias 


We require a homogeneous static magnetic field at the position of the magnetic 
sample in order to tune the frequency of our magnetic resonance. Tuning these 
resonances in the GHz regime requires fields on the order of several hundred 
millitesla at the mixing chamber stage of the cryostat at millikelvin temperatures. 


We utilize a pure-iron magnetic yoke construction to provide the field bias to 
our sample. The field itself is created by a homemade superconducting coil made 
out of a wire containing 54 niobium-titanium (NbTi) alloy filaments in a copper 
matrix with a total of about 4500 winding turns. Additional cylindrical samarium- 
cobalt permanent magnets, each yoke side equipped with one 10 mm x 5mm and 
three 10 mm x 1mm disks, create an offset field of about 178 mT, see Fig. 5.1(b). 
The magnetic field is concentrated around the magnetic sample inside the cavity 
resonator. We can tune the magnetic field with a rate of about 4.35mT/A, the 
actual current-field conversion factor changes between different cool-downs due 
to hysteresis effects in the permanent magnet and are also temperature dependent 
during a single cool-down. 
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20mm 


Figure 5.2: Clamping connection We use a clamping mechanism to create a low resistance, non-soldered 
connection between two superconducting wires. The inner copper clamps are electrically isolated from 
the outer holder by insulating sapphire wafers that are thermally conductive at low temperatures and 
pass the heat dissipated by the contact resistance to the cryostat plates. Insulating polyether ether 
ketone (PEEK) screws connect both parts. This mechanism is used at both pulse tube cooled stages 
where the cooling power is large enough to cope with the additional thermal input. (a) Closeup of the 
clamping mechanism. (b) Installation below the PT1 stage, connecting normal conducting copper wires 
to superconducting YBCO strips. 


With these tuning rates we require current values in the range of single digit amperes 
to perform our experiments. We use normal conducting 41.5 mm copper wires from 
room temperature to the PT1 stage that are connected to a stainles steel strip coated 
with yttrium barium copper oxide (YBCO), a high-T, superconductor provided to us 
by the group of Prof. Bernhard Holzapfel of the Institute for Technical Physics (ITP) 
at KIT. Since we cannot solder YBCO to the copper we utilize a clamping mechanism 
to press both wires together, leaving only a small contact resistance. The dissipated 
heat has to be transferred to the thermal anchoring point at the plate while at the 
same time maintaining galvanic isolation between current leads and cryostat ground. 
We use two interleaved clamps made out of copper with a small gap in between that 
is filled with several sapphire wafer pieces and insulating polyether ether ketone 
(PEEK) screws that connect both parts. Sapphire is an insulating material regarding 
charge transport but a well conducting material at low temperatures regarding 
thermal excitation. The clamps were designed and fabricated by our mechanical 
workshop and are depicted in Fig. 5.2. The same technique is employed at the PT2 
stage, connecting the YBCO wire to the low-T. superconductor NbTi wire with 
Te ~ 10K that is used for the magnetic coil. 
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5.4 Cryogenic setup 


Any further thermal anchoring at the still plate and intermediate plate is done by 
winding the wire around a copper post and gluing it with epoxy, with the insulation 
coating of the wire preventing an electrical short to ground. Between the lower 
plates, thermal conduction of phonons within the copper matrix is prevented by 
wet etching of the copper between the still and intermediate stage, as well as the 
intermediate stage and the mixing chamber. We etch the copper with with a 30 % 
nitric acid (HNO3) solution over a distance of several centimeters and glue the 
bare NbTi filaments with GE low temperature varnish to account for the reduced 
mechanical stability. This setup is tested up to currents of 10 A with the temperatures 
of the cryostat at the normal operating levels. 


5.4.2 Microwave setup 


In order to minimize the thermal load to the stages, we use coaxial cables of different 
material between the different temperature stages for the microwave wiring. We 
further minimize additional thermal noise input by attenuators and maximize our 
signal strength on the output side by amplifies and filtering. 


On the input side we use stainless steel cables all the way down to the base plate 
with thermal anchoring at all plates. All microwave components influencing the 
microwave signal such as attenuators, circulators, and amplifiers must therefore 
operate in this frequency range as well as the temperature of their local environment. 
Microwave attenuators are inserted at different temperature stages to reduce the 
microwave signal power and to avoid thermal excitation at higher temperatures 
plates. We insert 20 dB attenuators and thermally anchor them at the PT2 stage, 
still stage, and base plate at typical temperatures of Tpr2 = 3K, Tsin = 730mK, 
and Tpase = 25 mK to reduce the Johnson-Nyquist noise power by a factor of 100 at 
each stage. We attribute additional 15 dB attenuation to the microwave cables and 
connections at our measurement frequencies prior to the sample input. 


The reflected measurement signal from the sample is separated from the input 
signal by two cryogenic circulators connected in series. We therefore avoid the signal 
carrying the measurement data to be attenuated by the incoming lines and are able 
to further amplify the signal. On the outgoing side, we avoid cable losses by using 
superconducting niobium coaxial cables from the mixing plate to the PT2 stage. 
Additionally, the circulators shield the sample from unwanted incoming electric 
noise created by the amplifiers. We use two Quinstar QCY circulators operating 
between 4 GHz and 8 GHz that are shielded against 100 mT of magnetic field. At the 
PT2 stage we band-pass filter the signal by a MiniCircuit VHF-3100+ limiting the 
measurement to a band between 3.4 GHz and 9.9 GHz, before it is amplified by a Low 
Noise Factory LNF-LNC4_8C s/n 301C high electron mobility transistor (HEMT) 
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amplifier. The HEMT has a dedicated amplification band between 4 GHz and 8 GHz 
with an amplification of 40 dB and an added noise equivalent temperature of 2.0K. 
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Figure 5.3: Schematic setup for spectroscopic and time resolved measurements Depending on the 
type of experiment we either use the spectroscopy setup or the time domain setup (red) to measure the 
sample in the cryostat (blue). The components in between sit at room temperature and are used as part 
of both measurement schemes. 
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5.5 Spectroscopic measurements 


We use copper-nickel wires for the connection from the PT2 stage to room temperature. 
A second amplifier at room temperature is used to amplify the signal again 
before being recorded. We use a broadband Narda-MITEQ AMF-5D-00101200- 
23-10P amplifier with a gain of 40dB in an operating frequency range between 
1 GHz and 12 GHz. 


5.5 Spectroscopic measurements 


In spectroscopic measurements we probe the sample with continuous microwave 
signals from a Keysight E5071C vector network analyzer (VNA) that measures the 
complex scattering matrix element S21. While we measure the reflection characteristics 
of our sample, the wiring through the cryostat is built to appear as a transmission 
signal for the VNA. By analyzing the amplitude and phase response we extract 
resonator parameters or measure the hybrid systems by varying the excitation power 
or the applied magnetic field. An additional microwave tone, for example to measure 
the AC stark shift of the qubit, is provided by an Anritsu MG37020A microwave 
source and is fed into the input line by a directional coupler. The spectroscopic 
measurement setup is schematically depicted in Fig. 5.3. 


5.6 Time resolved measurements 


The time-dependent properties of the qubit are measured by manipulating its 
state with pre-defined microwave pulses and record the response of the system 
after a free evolution time At. The pulses rotate the Bloch vector around different 
axis of the Bloch sphere, while the final measurement always is an averaged qubit 
state measurement. 


We use a field programmable gate array (FPGA)-based setup developed by the group 
of Prof. Marc Weber of the Institute for Data Processing and Electronics (IPE) at 
KIT [184]. The FPGA board handles the pulse sequence and trigger as well as the 
data recording, while an analog microwave front-end is used to ensure synchronized 
pulses, amplification, and attenuation as well as frequency conversion to and from 
gigahertz frequencies. A schematic sketch of the setup is depicted in Fig. 5.3. 


The FPGA hardware platform runs a Linux operating system that eliminates the 
need for an external control computer during the experiment. With the whole 
experiment, data collection, and data processing being handled by the FPGA we are 
able to shorten the overall experimental time compared to traditional setup using 
an arbitrary waveform generator. A sequencer controls the pulse generation and 
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triggering for manipulation, readout, and data recording. It operates at 125 MHz clock 
frequency, leading to a time resolution of 8 ns. For convenience, the pulse sequences 
of frequently used qubit characterization experiments are pre-programmed, with 
only the time delays as input parameter. The sequencer triggers the pulse generator 
that feeds the amplitude modulated pulse to a digital-to-analog-converter (DAC) as 
the intermediate frequency (IF) for further analog signal processing. The MHz IF 
signal is mixed by an I/Q mixer with a local oscillator (LO) GHz signal, resulting 
in the sum and the difference of both frequencies. This converts the signal to the 
qubit or readout resonator into the GHz regime. The response from the system 
with respect to the readout pulse is down-converted using the same power divided 
LO signal as the up-converted readout pulse. The signal is low-pass filtered, 28 dB 
amplified, digitized by the analog-to-digital-converted (ADC), and recorded by the 
FPGA module. 


We use a Xilinx Zync UltraScale+ MPSoC ZCU 102 evaluation board with a 
Zync UltraScale+ XCZU9EG-2FFVB11561 multiprocessor system-on-chip together 
with two Texas Instrument evaluation cards as ADC (ADS54J69EVM) and DAC 
(DAC39J84EVM). Two Anritsu MG37022A microwave sources are used as local 
oscillators for manipulation and readout. We use a Marki IQ-0318L mixer for the 
manipulation signal and a Marki MLIQ-0416L mixer for the readout and recording 
signal. Manipulation and readout input into the fridge are combined by a Marki 
C20-0R520 directional coupler where the readout is —20 dB attenuated and the down- 
converted recording signals are amplified by Mini-Circuits ZFL-1000H amplifiers. 
The LO signal for readout and recording is split by a Mini-Circuits ZFRSC-1235+ 
power divider. 


5.7 Calibration routine before fitting the 
reflection spectrum 


The measurement setup, with its cable damping, cable connections, and all other 
added microwave components, contributes to the overall measured signal of the VNA 
due to frequency dependent attenuation and combined standing waves because of 
impedance mismatches and reflections in the coaxial cables. In order to separate the 
signal from the setup and from the sample we estimate the magnetic field independent 
baseline signal and consider it during the fitting process of the reflection spectrum 
of the coupled system. The procedure is detailed for the magnon-photon hybrid 
system, but is also used for the tripartite system considering the three hybridized 
modes. At the same time we reduce the number of free fit parameters, where we 


50 


5.7 Calibration routine before fitting the reflection spectrum 


obtain numerical values from different fit procedures and fix these values in the fit 
of the adapted spectrum. 


5.7.1 Baseline estimation 


Hybridizing magnon and photon results in two resonance dips in the in the reflection 
spectrum. While the measured signal outside of the dips stems from the measurement 
setup, the dip signal is a combination of setup and system. We estimate the setup 
baseline by via the field dependent hybridized modes, with which we identify the 
probe frequencies which contain only background information for each magnetic 
field value. We track the resonance minima of both cavity magnon polariton branches 
and cut a frequency band with a width of three times the loaded resonator linewidth 
below and above the tracked dressed frequencies from the original data. These 
system responses are influenced by the underlying physical system. Due to the strong 
coupling in our system, the two polariton branches separate in the anticrossing 
enough from each other to allow access to the measurement background of the setup 
even at the original resonator frequency. We determine the background by mediating 
the respective data for each probe frequency wp over all occurring current values 
and create an estimated response signal of the measurement setup. 


5.7.2 Adapting the fit routine 


The frequency dependent background |S11|yack (wp) is included into the fitting 
process of the reflection spectrum to only fit the signal of the underlying physical 
system. This is realized either by normalizing the original data and fitting it to 
Eq. (2.59) or by modifying the fit function and fit it to the original data. We adapt the 
fit function for the amplitude of the reflection spectrum to 


[S11) fp (cop) = [S11] (wr) - [Sti lpack (WP) 


14 Ke 


‘|Sitlback (WP). 6-1) 


i (Ww; — wp) + 3+; 


Since |S11|pack (Wp) is a single, known number for each wp we do not increase the 
number of fit parameters. 


We further reduce the number of fit parameters from six to only one, the internal 
magnetic linewidth km, by determining system parameters from different measure- 
ments or data analysis. The bare cavity resonator frequency wr, the field dependent 
magnon frequency Wm (I) and the coupling strength g are determined by fitting the 
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tracked minima of the two branches to the eigenvalues w of the 2 x 2 coupling 


matrix in the set of equations describing the two coupled harmonics oscillators, 
Eq. (2.40). The two resonator linewidths x, and x. are obtained by the zero-current 
measurements of the cavity resonator, table 6.1, which are independent of the applied 
magnetic field. 


5.8 Measurement software: qkit 


We perform our experiments with qkit, an open-source Python software suite 
developed at KIT [185]. The main focus lies on standardized measurement routines 
such as sweeps of multiple parameters in spectroscopic measurements and an 
integration of time-resolved qubit measurements. Instrument drivers handle the 
device communication and make the measurement code independent of the specific 
devices. This makes measurements flexible and adaptable for new users and new 
measurement ideas. 


The software has its roots, in particular the basic device communication, in the 
QTLab software package’, written by R. Heeres. We adapt his approach to our needs 
and the software is publicly available in a repository on the GitHub development 
platform. Our main focus is to build a common platform for all measurements in our 
lab with a powerful dynamical data viewer based on the PyQtGraph package” and 
an efficient data storage using the HDF format*. Qkit now provides a measurement 
environment for microwave and transport measurements, typically controlled via 
Jupyter Notebook?. It features data acquisition, storage, and visualization as well as 
some basic analysis tools for microwave resonators. 


2 QTLab, http: //www.pyqtgraph.org - accessed Oct. 18, 2019 

3 PyQtGraph, https: //www.pyqtgraph.org - accessed Oct. 18, 2019 
4 HDF Group, https: //www.hdfgroup.org - accessed Oct. 18, 2019 
5 Project Jupyter, https: //jupyter.org - accessed Oct. 18, 2019 
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6 Cavity resonator characterization 


In this chapter we present the microwave spectroscopy characterization of the two 
cavities used throughout this thesis. This resonance characteristics are relevant not 
only for fitting parameters, but also define the number of photon excitation of the 
system. For both cavities we focus especially on the resonance spectrum of the TE102 
mode, since the photons exciting these resonances are used as a component in our 
hybrid systems. The additional TEıoı and TE103 modes are only used for second 
tone excitation or readout purposes, where only the resonance frequency is relevant. 
We take a closer look to the mk_3 cavity resonances, used in experiments without 
qubit, with respect to the excitation power. This allows us to attribute any power 
dependent changes to the corresponding component of the hybrid system. 


6.1 mk_3 cavity charaterization 


We characterize the TE;92 resonance of the mk_3 cavity at T = 55 mK for various 
probe powers corresponding to the experiments in Sec. 7.2 by using a circle fit 
routine. The cavity cutout responsible for the resonance has a spacial volume of 
83 x 40 x 3mm?. 


Table 6.1: Resonance parameters cavity mk_3 TEıo2 resonance mode. The table includes the zero-current 
resonator frequency wi? and the loaded, coupling, and internal quality factors Qı, Qe, and Qi. 


wr °/2n (GHz) Qi Qc Qi 
5.239474 + (2-1076) 3084.445.1 5438.9+41 7124.9 + 34.3 


Since the magnetic part of the measurement setup features permanent magnets, we 
already apply a magnetic field to the YIG sample at any time that slightly dresses 
the bare resonator frequency wr. At zero applied current we measure a dressed 
resonance frequency wi?) With a probe power P = —140dBm, the fit determines a 


zero-current resonance frequency w!~° = 5.239 474 GHz + 2 kHz and quality factors 


53 


6 Cavity resonator characterization 


Qı = 3084.4 + 5.1 MHz, Q. = 5438.9 + 4.1 MHz, and Q; = 7124.9 + 34.3 MHz, see 
Figs. 6.2(a), (c), and (e). We compare the measured resonance frequency with finite 
integration simulation using CST Studio Suite 1 which results in a simulated 
wsim — 5,212 GHz that matches well the measurement. Besides an initial change of 
the quality factors we observe the TE102 resonance parameters to be constant with 
the applied power, see Fig 6.1. 


a 3098 — - - r z d - r - - r 
= 31 4 
© 3094} 4 
S N 2} ] 
3090} 1 3 
£ = ot | 
= Š 
3 3086 3 ol | 
g < 
3 3082+ d i 
3078 — f , i f j ‘ , ‘ y 
b 54457 r 7 r r 
© 5440} | ] | 
Re 
o 
o 
E 5435} ] è se Fd 
= 
D 
2 54301 i J 8 Y 
= 5 Y 
a N S 10 y4 
3 54251 4 } + + y € j 
t S ¿01 F 
5 10*t y 
5420 : Ea y 
7240 ~ 1 r r 7 x a 
c o 103$} Y 
= T je 
[67 2 ‘ 
5 7200F 1 gl # 
g = i 
£ e Ps 
© 7160+ 4 2 z 
» z 10; Y 
© Fg 
5 # 
£ 7120} 4 it # 
Y 
08020 -120 -100 -80 -60 -140 -120 -100 -80 -60 
Probe power P (dBm) Probe power P (dBm) 


Figure 6.1: Resonator parameters of the TEıo2 resonance of the mk_3 cavity (a-c) Loaded, coupling, 
and internal quality factors of the cavity resonance against probe power. The data was taken at 
T = 55mK with zero current applied to the magnetic coil and does not show a power dependent 
behavior. (d) Shift of the fitted cavity frequencies Aw = [wr (P) — wr (P 140 dBm)] with compared 
to the measurement at lowest probe power at zero current. Similar as with the quality factors, the cavity 
frequency does not show a power dependence. (e) Calculated average photon number in cavity against 
probe power. The fit shows a linear dependence of the photon number calculated with Eq. (2.15) to the 
input power. Note that this plot features a log-log scale, making the fit linear again. The errors on the 
average photon number are estimated to be smaller than 0.35% and are not visible in this plot. 


1 Dassault Systèmes Simulia, https: //www.3ds.com - accessed Oct. 18, 2019 
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6.2 mk_4 cavity charaterization 


We estimate the average excitation (n) in the cavity corresponding to the applied 
power with Eq. (2.15) and obtain 


(n) = 62.046 - P}.0008 gwt. (6.1) 


The single photon parameters of the TE102 resonance of the mk_3 cavity at T = 25 mK 
are also listed in table 6.1. In addition, we obtain the TE101 and TEj93 resonance fre- 
quencies at (0; 39, = 4.174811 GHz + 5 kHz and w! 703 = 6.644516 GHz + 8 kHz, re- 
spectively. From simulation we obtain the resonance frequencies watt = 4.162 GHz 


and am. = 6.606 GHz. All simulations match the measurements with an error of 
less than 0.6 %, with the simulated values all being larger than the measurements. 


6.2 mk_4 cavity charaterization 


The mk_4 cavity is used for the experiments on the strongly coupled magnon- 
photon-qubit system. Since the resonance frequency of the TE102 almost matches the 
first qubit transition frequency, both systems hybridize and form dressed state at 
low powers. 


Table 6.2: Resonance parameters cavity mk_4 TE102 resonance mode. The table includes the zero-current 
resonator frequency wi=0 and the loaded, coupling, and internal quality factors Qı, Qe, and Qi. 


wr ?/2n (GHz) Qı Qc Qi 
5.453271 + (8 : 1076) 1639.4 +3.0 2076.4+2.2 7790.2 + 56.4 


The difference in resonance frequency compared to the mk_3 cavity is realized by 
changing the resonance volume to 83 x 37 x 3m3. We access the pure cavity resonance 
parameters by applying a large probe power that decouples both components. The 


fit determines a zero-current resonance frequency wł™? = 5.453 271 GHz + 8 kHz 
and quality factors Qi = 1639.4 + 3.0 MHz, Q. = 2076.4 + 2.2 MHz, and Qi = 
7790.2 + 56.4MHz with P = —80 dBm, Figs. 6.2(b), (d), and (f). We determine 
the TE101 and TE;93 resonance frequencies at wio = 4.468 080 1 GHz 0.8 kHz 
and cod = 6.799 4644 GHz + 0.4 kHz, respectively. From simulation we obtain 
the resonance frequencies wh = 4.337 GHz, witty = 5.345 GHz, and with = 
6.690 GHz. All simulations match the measurements with an error of less than 3 %, 
with the simulated values here all being smaller than the measurements. Looking 
at the resonance parameters of the mk_3 cavity, we assume constant resonance 
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parameters in this cavity with respect to power as well. The resonance parameters of 
the the TE102 resonance of the mk_4 cavity at a probe power P = —80dBm and a 
temperature of T = 25 mK are also listed in table 6.2. 
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Figure 6.2: Cavity resonance measurements We plot the reflection spectra of the TE102 resonances and 
plot the amplitude |Sı1| and phase arg (S11) against the probe frequency, and the spectra’s imaginary 
part Im [S11] against the real part Re [S11] for both the mk_3 ((a), (c), and (e) at T = 55mK and 
P = —140 dBm) and the mk_4 cavity ((b), (d), and (f) at T = 25mK and P = —80 dBm). The baseline 
difference of about 20 dB between plots (a) and (b) is the result of an additional attenuator in the system. 
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7 Magneticlinewidth measurements 


This chapter focuses on experiments on a coupled magnon-photon system as 
described in Ch. 2. A coupling of these two excitations results in a harmonic 
hybrid system. While we are able to decrease excess excitation by the temperature 
of the environment and can control the excitation power we use to probe the 
system, the excitation is still Boltzmann-distributed over all possible evenly spaced 
excitation levels. 


We focus on the internal magnon linewidth that we extract from the microwave 
reflection spectrum of the coupled system and sweep temperature and excitation 
power to identify the dominant loss mechanism in the different environment 
regimes. This is an important step towards more sophisticated utilization of 
magnetic components in combination with superconducting qubits, which are 
operated at millikelvin temperatures in the single photon regime. 


7.1 Linewidth from 300 K to 2.4 K 


In a cooperation with the University of Mainz, we first characterize our YIG sample 
in a different measurement setup at higher temperatures. We measure the magnetic 
linewidth of a YIG sphere from room temperature down to about 2.4K [76]. The 
measurements are carried out similarly to the ones at millikelvin temperatures 
described in more detail in Ch. 5, but at fixed microwave excitation power. 


The system consits of a YIG sphere with a diameter of d = 0.5 mm mounted in a 
reflection cavity with a frequency of 6.53 GHz. With a *He continuous flow cryostat 
we are able to sweep the temperature from 2.4 K to 290K in a single measurement 
run and apply the required static magnetic field of 227 mT with a superconducting 
coil. The exact experimental details and measurement data can be found in Ref. [76]. 
The evaluated data for the magnon linewidth km is shown in Fig. 7.1. 


We observe an increase of the magnon linewidth with increasing temperatures 
from km /2rt = 1.16 + 0.14 MHz at T = 2.4K towards a dominant peak at about 40 K 
and decreasing with a relatively broad second maximum around 100K to about 
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Km/2m = 1.61 + 0.25 MHz at room temperatures. The main contribution to the 
linewidth in this temperature range is attributed to slow-relaxing rare-earth ion 
impurity scattering [119, 123, 124], see Sec. 2.5.1. 
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Figure 7.1: Magnon linewidth for temperatures from 2.4 K to 290 K We model the linewidth with two 
separate curves considering two rare-earth impurity scatterers with temperatures AQ, /kg = 43K and 
hQz/kg = 90K (dashed red lines). These contributions sum up to the blue line. The measured magnon 
linewidth at single kelvin temperatures is attributed to two-magnon scattering processes at the rough 


sample surface. The estimated contribution of Kgurface/2™ ~ 1.1 MHz compares well to the measured 
1.61 + 0.25 MHz at 2.4K. 


Using Eq. (2.61), we are able to model the linewidth data with two separate RE scatters 
with temperatures of AQ1/kg = 43K, a linewidth of 7 MHz and AQ2/kg = 90K 
with a linewidth of 3.5 MHz (red dashed lines), and combine them to the blue line, 
which matches our linewidth data well. 


At the lowest temperatures our data deviate from the curve modeling the slow- 
relaxing impurities that vanishes at 10 K. We attribute this to contributions from 
two-magnon scattering processes at the rough sample surface [125-127], see Sec. 2.5.2. 
We use Eq. (2.63) to model the influence of the surface scattering to our linewidth 
data. We describe the surface by hemispherical pits with a size of 2/3 of the 
roughness of the polishing material [128], rpit = 2/3 x 50nm. We determine the 
saturation magnetization to be 4yMs = 282 mT. This results in a surface scattering 
contribution Of Ksurface/ 27 ~ 1.1 MHz, matching well our measured linewidth of 
Km / 2m (T = 2.4K) = 1.61 MHz. 
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7.2 Cavity magnon polariton in the quantum regime 


Coil current I (A) 
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Figure 7.2: Cavity magnon polariton in the quantum regime of excitation (a) We plot the absolute 
value of the reflection spectrum |S1| (color coded) with respect to the probe frequency and the applied 
current at a temperature T = 55mK and a probe power P = —140dBm. The resonance dips show the 
dressed photon-magnon states, which form an avoided level crossing towards a frequency degeneracy 
point at Io = 2.09 A, corresponding to an applied magnetic field of By = 186.98 mT (dashed vertical 
line). The inset displays the squared gradient of the zoomed-in amplitude data. The kink in the data 
represents a weakly coupled magnetostatic mode. (b) Raw data of the cross section at the center of the 
avoided level crossing (black) and fit to input-output formalism following an adapted fit procedure 
(green). The data is normalized by the field independent background before fitting and the background 
is multiplied back to the fit to display it over the raw data. The fit yields a magnon linewidth of 
Km/2rt = 1.82 + 0.18 MHz at these conditions. 


7.2 Cavity magnon polariton in the quantum regime 


We use a (different) YIG sphere with a diameter d = 0.5mm together with the mk_3 
microwave cavity for further investigation on the loss mechanism of magnons in 
the quantum regime of excitations. The color-coded amplitude signal of a typical 
microwave reflection spectrum of the hybridized magnon-photon system for different 
magnetic field values is displayed in Fig. 7.2(a). The data was taken at a temperature 
of 55 mK with an excitation power of —140 dBm at the SMA input port of the cavity 
resonator. This power corresponds to an average number of (n) = 0.6 total excitations 
in the system, shared between magnon and photon. We observe a large anticrossing 
as the Kittel mode hybridizes with the cavity photon at constant frequency and we 
see the two branches of the dressed system. Due to the strong coupling with a mode 
splitting of 2g, we observe the two modes well outside of the frequency degeneracy 
point. In addition to the large anticrossing we observe a second smaller feature at 
slightly higher field strength, indicating the coupling of the photon mode to an 
additional magnetostatic mode within the YIG shpere. We enhanced the visibility of 
this feature by plotting the squared gradient of the data in the zoom-in insert plot. 
Since the second mode is only weakly coupled it does not form an anticrossing, 
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but only appears as a kink in the data. Such a magnetostatic mode is excited by 
a non-uniform ac magnetic field of the cavity resonance and are also reported in 
comparable experimental systems [61, 76]. 


Fig. 7.2(b) shows the raw data of the amplitude of the reflection spectrum (black) and 
the corresponding fit derived by input-output theory (green) for matching magnon 
and photon frequencies wm (Bo) = wr at an applied coil current of In = 2.086 A, 
corresponding to an applied external magnetic field of B (Io) = By = 186.98 mT. The 
current-to-field conversion derives from the known FMR dispersion relation for the 
Kittel mode in YIG and differs for different measurement dates and temperatures 
due to hysteresis effects in the permanent magnets. 


7.2.1 System parameters 


For the measurement depicted in Fig 7.2 we fit a bare resonator frequency 
w./2r = 5.23902 +2 - 107? GHz, a 452 kHz decrease compared to the zero-current 
measurement, see table 6.1. The magnetic field created by the permanent magnets 
slightly dresses the resonator frequency even at zero current, which leads to an 
increase compared to the bare frequency. The magnon frequency is given by 


wm (I) /2n = wi /2n + em -1 
kHz 


= 4.9817 + 2 - 1074 GHz + 122.0 +0.1 E 


È (7.1) 


which indicates a zero-current magnetic field of about Boff = 177.8 mT. We find 
a coupling strength g/2m = 10.39 + 0.17 MHz that is in good agreement with the 
theoretically predicted value of g™/27 = 12.48 MHz from Eq. (2.32). 


Table 7.1: Geometric system parameters For the theoretically predicted value of the coupling strength 
between magnon and cavity photon gi, we simulate different parameters with finite element simulations. 
The table includes the mode volume Vinode, the overlap factor y, the sphere diameter d, the corresponding 
number of participating spins Ns, and the spin number s. 


Vmode (m?) 7 d (mm) Ns 5 
5.406 -1076 0.536 0.5 1.37 1018 5/2 


We obtain the overlap factor y = 0.535, the ratio of mode volume in the cavity 
volume and the sample volume, used in the theoretical coupling formula from finite 
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element simulations using ANSYS HFSS !. For a sphere diameter d = 0.5 mm we 
assume a total number of participating spins Ns = 1.37 - 10!8 with a spin number 
s = 5/2. The geometric system parameters are also listed in table 7.1. 


We estimate the background signal strength from the measurement setup (Sec. 5.7) 
and fit the reflection spectrum at the degeneracy point, i.e. at matching magnon and 
photon frequencies to the modified formula derived by the input-output formalism, 
Eq. (5.1) and extract an intrinsic magnon linewidth km/2m = 1.82 + 0.18 MHz. The 
coupling strength exceeds both the total resonator linewidth x, and the internal 
magnon linewidth of km, thus being well in the strong coupling regime. While the 


linewidth is comparable to the value for room temperature measurements, it is 
significantly larger than the one at T = 2.4K [76]. The fit results for all relevant 
system parameters are also listed in table 7.2. 


Table 7.2: Cavity magnon polariton fit parameters The hybrid system parameters are determined by 
multiple measurements and fit routines to reduce the total number of free fit parameters per routine. 
The table includes the bare resonator frequency «wr, the zero-current magnon frequency ae", the 
coupling strength g, and the offset magnetic field created by the permanent magnets Boff- 


w,/2n (GHz) wi °/2n (GHz) g/2n (MHz) Bog (mT) 
5.23902 42-10? 4.9817+2:10°* 10.39 +0.17 177.8 


7.2.2 TLS losses in magnets 


In order to identify the dominant magnetic loss mechanism at mK temperatures 
we study the magnetic linewidth with respect to small temperature and power 
changes compared to the quantum excitation limit. At the degeneracy point we 
evaluate km at two power values, —140 dBm and —65 dBm in a temperature range 
between 55 mK and 1.8K and for detailed power sweeps with the input varying 
between —140 dBm and —65 dBm at two constant temperatures, 55 mK and 200 mK. 
The evaluated data are shown in Fig. 7.3. The cavity parameters and the coupling 
strength are not influenced by these temperature changes. By performing the data 
evaluation for each temperature stage on its own we compensate the hysteresis 
in magnetic field performance in the temperature sweep. The applied excitation 
power result on average in magnon numbers from (m) = 0.31 to (m) = 9.85 - 10°, 
with the corresponding scale also plotted in the the upper x-axis of Fig. 7.3(b). The 


1 ANSYS, Inc., https: //www.ansys.com - accessed Oct. 18, 2019 
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number of magnons is still much smaller than the number of participating spins 
on the order of 1018, therefore we expect a fully linear behavior of the magnons 
in contrast to intrinsic non-linear magnon behavior that is observed at magnon 
numbers comparable to the number of participating spins [58]. 


Temperature kg 7/h (GHz) Average magnon number (m) 
1 10 1 10 102 10° 10° 10% 10% 10 
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Figure 7.3: Magnon linewidth in the quantum regime of excitation (a) Temperature dependence of 
the magnon linewidth km for different probe powers. At low powers, the linewidth decreases with 
increasing temperature. We fit the tanh (1/T) trend from the standard tunneling model to the data, 
indicating that the linewidth is TLS loss limited in the quantum regime of excitation. The data at high 
powers (circles) do not follow this trend. The TLSs are already saturated by thermal excitation and the 
linewidth remains constant. (b) Power dependence of the magnon linewidth km for T = 55mK and 
200 mK. The temperature curves show a similar behavior. The linewidth stays constant for low probe 
powers and we observe a drop in km at powers of about —90 dBm. We fit a 1/,/(1+ P/P:) trend of the 
TLS model to both data sets and retrieve a similar critical power as well as matching xy and Ko¢ for 
both temperatures. All linewidth data shown here are extracted from the fit at matching frequencies. 


The magnetic linewidth decreases going to higher temperatures at small excitation 
power, but stays constant for the measurements at high excitation power, with the 
strong excitation linewidth matching the high temperature data at small excitation, 
Fig. 7.3(a). In the power sweep, km decreases at both temperatures with increasing 
excitation power, Fig. 7.3(b). We attribute this behavior to excitation losses into 
a bath of two level systems, see Sec. 2.5.3. The TLSs become saturated either by 
thermal energy at higher temperatures or resonant power absorption by the strong 
excitation. At strong excitation in the temperature sweep, all TLSs are saturated and 
Km Stays constant. We expand the function describing TLS losses, Eq. (2.73), by an 
offset linewidth ko, which describes the magnetic linewidth without TLS influence, 


tanh (hw;/2kpT) 
J1+P/P. 


Km (T, P) = Ko + Koff- (7.2) 
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From the fit, we determine the low temperature limit of the linewidth describing 
the TLS bath within the sample xo/2rt = 1.05 + 0.15 MHz and the critical power 
Pe = —81.0 + 6.5 dBm at the SMA port. We obtain the offset linewidth Koff/ 27 = 
0.91 + 0.11 MHz as a lower boundary for the internal linewidth without any TLS 
contribution in the saturation limit. We neglect the power dependence 1/ v1 + P/P. 
in the fit to the temperature sweep data in Fig. 7.3(a) since we operate in the single 
excitation regime with P < P.. The fit parameters to the TLS models are also listed 
in table 7.3. 


The high power data in Fig. 7.3(a) increases slightly at temperatures of about 1 K. 
We attribute this to the increasing contribution of scattering losses at rare earth 
impurities within the sample, which we also observe in Ref. [76]. 


Table 7.3: Fit parameters to the TLS loss model The table includes the low temperature limit of the 
magnon linewidth xg, its lower boundary without TLS contribution Koff, and the critical excitation 
power Pe. 


Ko /2m (MHz) kofg/2r7 (MHz) Pe (dBm) 
1.05 + 0.15 0.91 + 0.11 —81.0 + 6.5 


The linewidth decrease in the temperature sweep shows the thermal excitation of 
the resonator, where we recover the temperature corresponding to the resonator 
frequency, hwr = kgT,, with T, = 251.4 mK at half the TLS-induced linewidth change. 
The corresponding frequency scale to express the temperature is plotted on the 
upper x-axis in Fig. 7.3(a). 


With our measurements and the fits in Fig. 7.3, we can confirm incoherent coupling 
in a bath of TLSs as the dominant loss mechanism of magnons in the quantum limit 
of excitation and extend temperature dependent measurements [61] by including 
the power dependent measurements and analysis as well [62]. A similar behavior is 
also reported in thin film YIG structures [64]. 


7.2.3 Linewidth outside of the anticrossing 


So far, only the magnetic linewidth at zero frequency detuning was considered, 
where the energy in the system is equally shared between magnon and photon. We 
now want to observe the power dependence of the magnon linewidth at constant 
input power and instead control the energy in the magnon via frequency detuning. 
At constant input power, the magnon excitation depends on the frequency detuning 
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because fewer magnons are excited at off-resonant applied fields. Figure 7.4 shows 
the magnetic linewidth with respect to the applied magnetic field at temperatures 
T = 55 mK (subplot (a)) and T = 200 mK (subplot (b)). The corresponding magnon 
frequency is displayed on the upper x-axes and dashed horizontal lines mark the 
cavity resonance frequency. The different excitation powers are color coded and the 
legend is valid for both plots. We only display the highest excitation powers that 
show the decrease in linewidth. The x- and y-axis scaling is identical for both plots. 


a Magnon frequency w,/2x (GHz) b Magnon frequency w,,/2 (GHz) 
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Figure 7.4: Magnon linewidth outside of the anticrossing We plot Km against the external magnetic 
field for different probe powers at T = 55mK (a) and T = 200mK (b). The depicted probe powers 
correspond to the ones at the linewidth decrease in Fig. 7.1(b). The linewidth shows a global minimum 
in the anticrossing (dashed line), here we observe the strongest TLS saturation effects since the number 
of excited magnons is at its maximum. We observe local minima at 187.25 mT external magnetic field 
that correspond to the coupling to an additional magnetostatic mode within the sample (see inset of 
Fig. 7.2(a)). The linewidth increases with higher frequency detuning and reaches the TLS-influenced 
values outside of the crossing. 


For both temperatures we find a minimum of the magnon linewidth at matching 
frequencies at an external magnetic field of about 187 mT and the linewidth increases 
with the frequency detuning A = wr — Wm to match the low power values. 


At highest powers we observe for both temperatures a second local minimum at 
about 187.25 mT external field. We attribute this to the coupling to an additional 
magnetostatic mode within the YIG sphere, due to a non-uniform microwave 
magnetic field at the sample. This magnetostatic mode is also faintly visible in the 
color plot of the anticrossing. The inset plot in Fig. 7.2(a) shows the squared gradient 
of the zoomed-in amplitude data, where the small kink shows the weak coupling 
to the second magnetic mode. This accidental excitation of a magnetic mode has 
also been reported in similar experimental data [61, 76]. Coupling to a second mode 
increases the number of excited magnons and saturates more TLSs, which results 
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in the local minimum in the linewidth data. We cannot explain the instantaneous 
increase in linewidth in the T = 55 mK data in the vicinity of the anticrossing. 


We further quantify the energy ratio by which the Kittel mode is excited. Two 
resonating system that hybridize mix by a given amount depending on the frequency 
detuning. Since we only have direct access to the cavity, we calculate the energy 
stored in both branches as seen from the cavity and than weigh the energy ratio 
in the with the corresponding magnon /photon mixing ratio in each branch. Each 
anticrossing is separeted in its upper (+) and lower (-) branch, where we individually 
fit the resonances in each branch at the tracked dressed resonance frequencies with 
the circle fit algorithm [82, 83]. From the fit we obtain the total energy in units of 
excitations numbers for each branch, 


EL SI up, (7.3) 


with corresponding energy ratios 


(7.4) 


The used quality factors are the ones determined by circle fits of both Lorentzian 
dips and not the ones obtained by the pure cavity measurements. 


Each branch has a distinct mixing ratio of magnon and photon share, depending on 
the frequency detuning. This is reflected in the components of the eigenvectors X+ 
of the coupling matrix A in Eq. (2.39), 


(7.5) 


x. _(%2)__1_ NOFA 
u x?) J20 VOA E 


with Q = y4? + 4g?. The magnon/photon mixing ratio is given by the squares of 
the corresponding components of the eigenvectors. Due to symmetry reasons the 
magnon share m in one branch is equal to the photon share p in the other branch at 
any given A. Since we monitor the resonance mode of the cavity, we always observe 
a photon ratio larger than (or equal to) the magnon ratio. This leaves us with 


m+ =p = eai (7.6) 


ratio 


The total energy ratio in the magnon Er; p 


from both branches, 


is the sum of the magnon/photon energy 


ae + |2¢(2)[" grato (7.7) 
‘ 2 2 f 1 2 3 A 
Ea > 1x! | pratio 2b xp! Efatio =7= Ea (7.8) 
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The energy ratios for the measurements at —65 dBm excitation power are displayed 
in Fig. 7.5 with the same x- and y-axes and scaling for both subplots and representing 
the same magnetic field section as in Fig. 7.4. As expected, we observe the maximum 
magnon excitation ratio at zero detuning with 50% of excitation and a decrease 
to 20 % excitation share at the borders of the subplots. The coupling to the second 
magnetic mode and the increase in magnon number is also visible as a local maximum 
at both temperatures. 
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Figure 7.5: Energy ratio stored in the magnetic sample and the cavity photon The magnon energy 
ratio shows the coupling to the second magnetostatic mode at 187.25 mT. The ratio drops from the 
maximum at A = 0 to 20% at the plot boundaries. The x-axes scaling is the same as in Fig 7.4. (a) 
55 mK, (b) 20 mK. 


7.2.4 Possible origin models for magnetic TLSs 


The microscopic origins of TLS losses in magnetic material (and for that matter, 
dielectric TLSs as well) is not fully understood yet and part of ongoing research. 
There are multiple possible explanations why TLS losses are observed, such as 
magnetic TLSs, surface spins or magnon-photon coupling. 


Magnetic TLSs were proposed with an analog behavior to electric dipolar coupled 
TLSs [186-189] and measured in spin glasses by thermal conductivity, susceptibility 
and magnetization measuements at low temperatures [190, 191]. Amorphous YIG 
shows spin glass behavior [192] and therefore it is plausible to observe these effects 
in our crystalline YIG sample as well. With already observed rare earth impurities in 
our sample [76], it is likely that it contains structural crystal defects as well. This 
is again analogous to materials with electric dipolar coupled TLSs, where TLSs 
appear largely in disorderd crystals but were also detected in smaller densities in 
single crystals [193]. 
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Another possible origin of the losses could be surface spins, since they were observed 
as an important loss mechanism in cCQED experiments [147, 148] showing TLS loss 
behavior. In such a case we would find a changing number of participating spins 
in our system. We evaluate the coupling strength g to find a power or temperature 
dependence and observe such a change in Ng, see Fig. 7.6. We find an increase in 
the coupling strength of about 1% at the saturation conditions for the TLSs. With 
g x yN; this translates to an increase in the number of participating spins on the 
order of 2%, possibly due to an increased participation of spins on the surface that 
do not contribute to the coupling strength at low excitation. 
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Figure 7.6: Power and temperature dependence of the coupling strength The temperature dependence 
(a) and power dependence (b) of the coupling strength evaluated at the same conditions as shown 
in Fig. 7.3. At constant input power, the coupling strengths shows no conclusive temperature trend 
with little changes. At constant temperatures it increases by an order of 1% at high excitation powers. 
From these increase in coupling strength we obtain an increase in participating spins on the order 
of 2%. Due to the small calculation error, most error bars are not visible in the plot. The numerical 
values for the coupling strength at matching conditions in both subplots do not necessarily match, 
since these measurements were carried out during different cool-downs which influences the exact 
geometric shape of the copper cavity. 


A loss mechanism by magnon-phonon coupling and subsequent phonon losses due 
to TLS coupling is neglected because these magnon losses are expected to be much 
smaller than the Gilbert damping [194] for Kittel mode magnons with wave vector 
k = Din YIG. 
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In this chapter we introduce the transmon qubit as a nonlinear element into the 
experimental system. Coupling the qubit to the magnetic excitations allows us 
to individually address single magnons in the system. We explore two different 
coupling regimes with two different microwave cavities. First, we create an indirectly 
coupled qubit-magnon system that is located in the dispersive regime and read out 
via standard qubit measurements. In a second step, we match the cavity frequency 
to the qubit frequency and create a strongly coupled system of all three components; 
magnon, photon, and qubit. 


First, the qubit is characterized in the cavity with spectroscopic and time resolved 
measurements without any applied current. We then measure the coupled system 
in the dispersive limit and focus lastly on the strong coupling regime, where we 
perform power and temperature sweeps complemented by corresponding analytic 
fits and simulations. 


8.1 Qubit characterization 


We first characterize the Josephson junction in the transmon qubit with resistance 
measurements at room temperature. From spectroscopic measurements at mil- 
likelvin temperatures we obtain the qubit transition frequencies and determine the 
characteristic qubit energies Ej and Ec. We calibrate the excitation number in the 
TE103 mode by the ac Stark shift [16], a power dependent shift of the qubit frequency, 
to later operate our system in the single photon regime, and measure the dissipative 
dynamics with pulsed time domain measurements. All these measurements are 
performed using the mk_3 cavity without any applied current. 


Fig. 8.1 shows a false-colored micrograph image of the used transmon qubit chip. 
The design is largely dominated by the pads of the shunt capacitor on both sides of 
the Josephson junction. The junction itself a size of approximately 100 x 95 nm? and 
cannot be resolved in the optical images even at a 100x optical zoom level and is 
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sketched schematically into the picture. The normal state resistance of the junction 
Rn = 7.51 kO, is obtained with a 2-point probe setup at room temperature. Using the 
Ambegaokar-Baratoff relation between junction resistance and the critical current 


I. [195] we find 


TA 1 
lo = X R 71.11 nA, (8.1) 


with the superconducting gap energy of Al, A = 3.4 : 1074 eV. This critical current 
translates to a Josephson energy En = 35.3 GHz - h, which we combine with the 
simulated charging energy EM = 277.8 MHz - h to result in a characteristic energy 
ratio E; /Ec = 127.1 and an expected qubit transition frequency w00 /2rı = 8.58 GHz 
based on the room temperature measurements. Comparing these values with the 
measured qubit frequency and energy ratio measured at millikelvin temperatures 


and discussed below, where we find a significant deviation from the expected value. 


Figure 8.1: Transmon qubit chip False colored micrograph of the transmon qubit. The capacitor pads 
dominate the geometrical design of the structure, with only small connecting Al strips leading to the 
Josephson junction. The junction itself with an estimated size of 10095 nm? is too small to be seen here 
even in the zoom-in and is only depicted schematically. 


8.1.1 Spectroscopic measurements 


We install the qubit chip into the cavity box, cool the system down to 25 mK, and 
probe it with continuous microwaves using the spectroscopy setup to measure the 
qubit transition frequency and the coupling strength between qubit and cavity 
resonator. The qubit transition frequencies are observed as dips in the dispersive shift 
measurements of the TE103 readout mode frequency at wı03/2rt = 6.64 GHz with 
respect to the qubit probe frequency and probe power, see Fig. 8.2(a). We observe 
multiple dips appearing and broadening at increasing probe powers, of which we 
identify the two dips at the highest frequencies with the energetically lowest qubit 
transitions. The fundamental |0) — |1) transition at wo; /2m = 5.4315 GHz and the 
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two photon |0) — |2) transition at 5wo2/2m = 5.3526 GHz. With Eq. (3.24) we find 
for the characteristic energies 


1 
fe? (Eo - 5a) = 167.8MHz - h, (8.2) 
Eo + Ec)* 
Ey = EutEe) _ 99.355 GHz -h, and (8.3) 
SEC 
EJ Ea = 1992) (8.4) 


This positions our qubit well in the transmon regime. We identify the third dip in 
the dispersive shift spectrum at 5.252 GHz as the |1) — |2) and the three photon 
|0) > |3) transitions, which are expected at the same frequency. Based on the 
calculated parameters we would expect these transition at w 2/2m = 5.265 GHz, 
indicated by the corresponding dashed line in Fig 8.2(a). 


a b 


En 

P 

u 
tin 
N 
N 


Probe frequency wp/2x (GHz) 
e 
un 
u“ 
N 
o 


Probe frequency wp/2x (GHz) 


u 
N 
u 


5.18 
-115 -110 -105 -160 -140 -120 -100 -80 
Qubit probe power P (dBm) Resonator probe power P (dBm) 


Figure 8.2: Spectroscopic measurements, qubit transitions and qubit-resonator coupling (a) The first 
transitions frequencies are obtained by evaluating the dips in the dispersive shift spectrum of the TE103 
readout mode that is 4103/27 = 1.14 GHz detuned from the qubit. The transition frequencies and the 
non-equidistant frequency spacing determine the characteristic qubit energies; E; = 23.355 GHz - h, 
Ec = 167.8 MHz - h, and Ej/Ec = 139.2. The calculated |1) — |2) transition matches well with the 
corresponding signal. (b) The TEj92 cavity mode (coupling mode) becomes dressed by the qubit 
presence (low power). At higher powers the qubit is saturated and the system decouples. The shift of 
nearly 20 MHz relates to a coupling between qubit and coupling mode of gq102/27 = 66.26 MHz. 


The measured fundamental transition frequency and the characteristic energies differ 
significantly from the room temperature measurements. The simulated charging 
energy is 65 % higher than the measured Ec and the Josephson energy obtained from 
the critical current measurement is 51 % higher than the value of the millikelvin 
measurement. With the relative deviations similar for both energies, the characteristic 
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energy ratios match within 10 % accuracy. This results in a fundamental frequency 
that is 3.15 GHz lower than initially expected. However, this does not play a role in 
the usability of the qubit chip for further experiments. 


We measure the coupling strength gq 102 between qubit and the TE102 coupling mode 
via its Lamb shift. The cavity resonator frequency is influenced by the presence of 
the qubit, similar to the Lamb shift in atoms. The qubit-cavity coupling gq i102 results 
from the frequency detuning between qubit and resonator mode, A192 = w102 — w01, 
and the frequency shift x between the coupled and uncoupled resonator. We probe 
the system with varying powers around the cavity resonance to saturate the qubit 
and measure the bare resonator frequency, see Fig. 8.2(b) and find the resonator mode 
power shifted by x192/2m = -19.9 MHz. With the frequency difference between 
cavity mode and qubit 4102/27 = -220.6 MHz we find [163] 


2 
x= A => 8410/21 = V/X1024102/27 = 66.26 MHz. (8.5) 


These measurements position our qubit in the dispersive regime with respect to the 
cavity resonances, where the frequency detuning between any resonator TE10p mode 
is much larger than its coupling strength with the qubit, A10p > 87,10p- This allows 
for the system to be described with the Jaynes-Cummings Hamiltonian, Eq. (3.29). 


8.1.2 Excitation number calibration 


In contrast to the estimation in the magnon-polariton case, Sec. 6.1, we are here 
able use the ac Stark shift, a photon number dependent shift of the qubit frequency, 
to exactly calibrate the number of excitations (n) in the TEj93 mode [16, 196-198], 
which is used as a reference for further experiments. 


The approximately diagonalized Jaynes-Cummings Hamiltonian, 
Be ae oa 1 
Om Ut /n~ 5 (w +x+ 2xa*a) 0 +â â, (8.6) 


2 
contains a mixed term &.0;4'4 that can be interpreted as a photon number dependent 


contribution to the qubit frequency, with 


(8.7) 


describing the ac Stark shift of the qubit frequency. 


We measure the qubit transition frequency versus the applied drive power at a probe 
power P = —140 dBm, see Fig 8.3(a). At an increased readout power we observe a 
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shift of the dressed transmon levels away from the readout mode, in our case to lower 
frequencies. We track the shifted qubit frequencies and fit the linear dependence on 


the excitation number, 
n103) = 12.69 - P22 4w-t (8.8) 
in 


with the coupling of the qubit to the TE;93 readout g103/2rt = 70.05 MHz and a 
frequency difference 4103/27 = 1.14GHz, see Fig. 8.3(b). This is a reduction of 
about a factor of 4.9 compared to the estimation of the population of the TE102 mode, 
Chap. 6. We attribute this to comparable Q factors at a higher resonance frequency, 
which leads to an overall smaller excitation number for the higher cavity mode. 
Based on the fit we perform the spectroscopic experiments in the dispersive limit at 
a readout power of Pp = —130 dBm corresponding to a single photon excitation. 
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Figure 8.3: AC Stark measurement (a) By increasing the drive power on the readout resonator we 
populate the cavity with more photons. The qubit frequency becomes dressed and shifts away from the 
readout frequency. (b) We track the shift in qubit frequency in the AC Stark measurement and calculate 
the exponentially power dependent photon occupation in the cavity. This improves the accuracy 
compared to the estimation done in the magnon polariton experiments. 


8.1.3 Time resolved measurements 


We obtain the dissipation dynamics of the qubit, the energy relaxation time T; and the 
decoherence time T7, by pulsed measurements using the FPGA-based measurement 
setup described in Sec. 5.6. 


We measure the energy relaxation time Ty by applying a previously calibrated 
n pulse with a length of Tun = 161.0ns that excites the qubit from |0) to |1) and 
measuring the qubit state after different waiting times At, see Fig. 8.4(a). Fitting an 
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exponential decay to the averaged qubit inversion, a measure for the probability to 
find the qubit in the excited state, we obtain Tı = 402.9ns. 


a b 
1} Tı = 402.9 ns y 11 TZ = 521.8ns 7 
= E 
¿2 2 
n Y 
v © 
2 E 0.5} 
2 2 
5 5 J 
o o readout 
inh AW 
ot L 1 1 1 A J ot 1 s L À 1 1 
0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 
Delay At (us) Delay At (us) 


Figure 8.4: Measurement of the qubit lifetimes (a) Applying a previously calibrated n pulse rotates 
the qubit into the excited state. We measure the qubit state after different waiting times Af and obtain 
an energy relaxation time Tı = 402.9 ns. (b) The decoherence time Tf = 521.8ns is measured with a 
Ramsey experiment, where the qubit state is rotated in a superposition state, evolves freely for times At 
on the equatorial plane, and is then rotated further around the x-axis. The applied pulse is frequency 
detuned from the qubit frequency by 5 MHz. 


We measure the dephasing time TF by a Ramsey pulse sequence, where we excite 
the qubit by a n /2 pulse in a superposition state on the equator of the Bloch sphere 
and let it evolve freely for times At without additional echo pulses, before rotating it 
by the same pulse and measuring the qubit state, see Fig. 8.4(b). 


Table 8.1: Transmon parameters The table includes the first qubit transition frequency wq, the 
anharmonicity a, the characteristic energy ratio E; /Ec, and the characteristic qubit times Tı and Ty. 


Wq/2m (GHz) a/2n (MHz) Ej/Ec T (s) Ty (s) 
5.4315 —179.5 1392 402.9 521.8 


The excitation pulse is detuned from the qubit frequency which results in Ramsey 
fringes [199], since the detuned tt /2 pulse adds an additional oscillation to the qubit 
state that decays exponentially. Fitting our data, we recover the 5MHz detuned 
excitation pulse and extract a dephasing time TY = 512.8 ns. All transmon lifetimes 
are also listed in table 8.1. 
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We estimate the losses due to spontaneous Purcell emission via the cavity linewidth 
Kr into the connected feed line to be Tp = 1/ yp = 1.138 ps with 


(8.9) 


While the 220 MHz detuning of the qubit to the TE102 coupling mode is quite small 
and leads to a relatively strong Purcell effect, our measured qubit lifetimes do not 
seem to be limited by these radiation losses. 


8.2 Magnon-vacuum Rabi splitting in the 
dispersive limit 


We create the quantum hybrid system between the Kittel mode and the supercon- 
ducting transmon qubit in the dispersive limit with the magnon resonance frequency 
matching the fundamental qubit transition. The nonlinearity of the qubit allows 
us to address the individual energy levels of the system and access single magnon 
excitation, in contrast to the previously discussed purely classical experiments con- 
necting magnon excitation with harmonic photon modes, where we simultaneously 
excited all energy levels with a Boltzmann distribution. While the direct coupling 
between the components is negligible, both couple to the TEj92 coupling mode, 
which leads to an effective coupling via virtual photon exchange between qubit 
and magnon. The system stays in the dispersive regime with both frequencies far 
detuned from the coupling mode. 


We tune the magnon frequency towards the qubit transition and monitor the 
dispersive shift of the TE103 readout mode while driving the system with a second 
tone in the vicinity of the qubit frequency in the single photon regime. The dispersive 
shift features a splitting in the qubit transition frequency that is attributed to a level 
hybridization due to a coherent exchange of excitation of the qubit and the Kittel 
mode with an effective coupling strength goa m /2m = 3.35 + 0.16 MHz, see Fig. 8.5. 
The coupling is mediated mainly by the TE102 mode, with an effective interaction 
Hamiltonian in the corresponding rotating frame, see Sec. 4.1, 


Ĥint = ng, (ato. + 90+) , with 


sth _ y £m10p8q,10p 
Sqm 3 W10p — Wg 


Due to the unwarranted effort to create the needed magnetic field strength to 
measure gm,103, applying a current of nearly 14 A, we only consider the contributions 
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of the TE101 and the TE102 mode to predict the effective coupling. With the measured 
coupling strengths connecting the individual components to the resonator modes 
we find 

¿05 mm ay 2 = -38.80 10.39: 66.26 


= 1265.1 ` 220.6 ) MHz = 3.30 MHz. (8.10) 


We do expect the approximation to be valid, given a TE101 contribution to Zqm 
of only 0.18 MHz = 5.4%, and the estimated coupling strength matches well the 
measurement. The expected effective coupling strength is found as half the distance 
of the hybridized modes in Fig. 8.5(a) but, given the participating linewidths 
kr /21t = 850 kHz, km/2r = 1.8 MHz, and xq/2n = T = 1.95 MHz, the effective 


coupling does not place the system in the strong coupling regime. 
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Figure 8.5: Magnon-vacuum Rabi splitting We hybridize the magnon and qubit excitation and measure 
the resulting magnon-vacuum induced Rabi splitting in the dispersive shift measurement of the driven 
system. The coupling strength is increased by exchanging the 0.5 mm YIG sphere (a) with a 1.0 mm 
YIG sphere (b), but the system does not stand in the strong coupling regime. 


In order to increase the effective coupling we install a larger YIG sphere witha diameter 
of 1.0mm, compared to the previously used 0.5 mm sphere, into the cavity. This 
results ina factor of 23 = 8 times more participating spins and an expected increase of 
the effective qubit-magnon coupling by a factor of v8 to a mm /2r = 9.34 MHz, 
which would position the system in the strong coupling regime. While the increased 
effective coupling strength of Y 2m = 8.21 + 0.11 MHz is well visible in the the 
dispersive shift measurement in Fig. 8.5(b) and matches the expected value with an 
error on the order of 10 %, the overall signal to noise ratio is too weak to unambiguously 
show two resonance minima for one magnetic field value. Nevertheless, we clearly 
observe the effective qubit-magnon coupling in the dispersive limit and are able 
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to confirm the square root dependence of the coupling strength on the number of 
participating spins. All coupling strengths are presented together in table 8.2. 


Table 8.2: Effective qubit-magnon coupling strength in the dispersive limit The table includes the 
measured effective qubit-magnon coupling strength ¿qm and the theoretically expected value En /2m 
for two different YIG sample. 


diameter (mm)  ¿qm/2rt (MHz) a /2r (MHz 
0.5 3.35 + 0.16 3.30 
1.0 8.21 + 0.11 9.34 


8.3 Strongly coupled tripartite system 


With the previously measured qubit-magnon hybrid system in the dispersive regime, 
the next step is to bring all three participating systems on resonance with each other 
to create a strongly coupled system consisting of magnon, qubit, and photon. 
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Figure 8.6: Strongly coupled qubit-photon system Power sweep of the reflection spectrum (a) and the 
dispersive shift measurement (b). The nearly frequency degenerate qubit-resonator system hybridizes 
and the dressed frequencies Wq + are observed at low power excitation. Increasing the power saturates 
the qubit and only the cavity resonance is observed. We obtain a qubit frequency wq/27 = 5.4304 GHz 
and a qubit-resonator coupling strength ggr/2rt = 67.65 GHz 


With the magnon representing a spin system, the qubit being an artificial atom, and 
the photon being a light quantum we therefore create a strongly coupled hybrid 
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system of light, matter, and spin. At the same time we increases the overall signal 
strengths and improve the signal to noise ratio by decreasing the frequency detuning 
between the photon coupling mode and the other two components. This is the 
next step towards more sophisticated implementations of magnetic components in 
quantum circuits that require coherent excitation exchange. 


We first characterize the qubit in the cavity resonator in the strong coupling regime 
with respect to its frequencies and coupling strength, then tune the Kittel mode 
frequency on resonance to create the tripartite system. Different probe powers allow 
for the qubit-resonator coupling to be turned on and off to individually study the 
resonator-magnon coupling. The system parameters are determined by the two 
component measurements and the full tripartite measurement, where we compare 
the obtained frequencies and coupling strengths. We analyze the spectrum with 
a three-component analytic function based on the input-output formalism and 
simulate the spectrum and the eigenenergies of the full tripartite Hamiltonian with 
qutip [200], a open-source python package for simulating the dynamics of open 
quantum system. We study the thermal influences on the system and compare the 
simulated spectra to the data measured at the corresponding temperatures. 


8.3.1 Strong qubit-resonator coupling 


Since we work with a non tunable qubit, we adapt the cavity resonator frequency by 
installing all components in the mk_4 cavity, whose TE102 coupling mode is designed 
to match the previously determined qubit transition frequency. The cavity resonance 
characteristics are discussed in Sec. 6.2. The cavity resonator hybridizes with the 
qubit and we observe hybridized qubit-resonator modes, in contrast to the Lamb 
shifted mode in the dispersive regime due to the presence of the qubit. Figure 8.6 
shows a power sweep of the reflection spectrum measured directly with the VNA (a) 
and the corresponding dispersive shift measurement of the TE103 readout mode of 
the driven system (b) without applied current. The qubitis populated with increasing 
power and decouples from the resonator. We observe the two hybridized modes in 
the low power regime and the bare resonator frequency is observed at high probe 
powers. In the intermediate regime the excitation is distributed into multiple modes 
that are averaged out in the spectroscopic measurement. 


The general scheme of the tripartite system is depicted in Fig. 4.1. We find 
two hybridized qubit-photon modes at Wgr,—/2nm = 5.373 GHz and war, + /2n = 
5.510 GHz in the single photon regime and observe the bare cavity resonance at 
wr /2T = 5.453 GHz. The dressed frequencies are described analog to Eq. (2.40). We 
approximate the qubit in the weak excitation limit to only have a single transition 
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and are able to determine the qubit frequency wo, = wg and the qubit-resonator 
coupling gar, 


Wq = War + + War, — Wr = 5.4305 GHz - 27 and (8.11) 


gar = for, + ww- — ww — w? = 67.60 MHz - 27. (8.12) 


The dispersive shift measurement in subfigure (b) confirms the measurement of 
the reflection spectrum obtained directly by the VNA, with additional features in 
the transition regime between coupled and uncoupled system that we attribute to 
dressed states with higher qubit modes. 


Coil current I (A) b Coil current I (A) 
3 4 5 6 3 4 5 6 


Probe frequency wp/2x (GHz) 


5.4 5.5 : : 5.4 5.5 
Magnon frequency wm/2r (GHz) Magnon frequency w,/2 (GHz) 


Figure 8.7: Current sweep of the tripartite system (a) We create the tripartite system with three 
resonance frequencies by tuning the Kittel mode frequency through the dressed qubit-photon states 
and fit the system by a constant cavity resonator frequency and two linearly changing frequencies for 
the Kittel mode and the qubit. The additional two anticrossing at the tripartite points have effective 
coupling strengths ¢/2m = 18.71 MHz and §/2m = 23.07 MHz.(b) Decoupling the qubit from the 
system by applying a probe power P = —95dBm recovers the magnon-polariton measurement similar 
to Sec. 7.2 and we obtain a resonator-magnon coupling strength of grm = 30.41 MHz. 


Comparing the system parameters with the previous measurements in the dispersive 
case using the mk_3 cavity, we find the qubit frequency to be shifted by —10 MHz. 
This is attributed to the temperature cycle of the sample during warm-up and cool- 
down, as well as the exposure to oxygen. The cavity coupling mode nearly matches 
the qubit frequency, being detuned only by A102 = 22.8 MHz. The qubit-resonator 
coupling remains nearly the same, with an increase of only 1.39 MHz despite a 
change in the effective qubit position with respect to the electric dipole of the cavity 
resonance field. 


As sketched in Fig. 4.1(c) we create the tripartite quantum system by tuning the 
Kittel mode into resonance with the two hybrid qubit-resonator modes in the 
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single photon regime at a probe power P = —145dBm. The measurement of the 
current sweep is depicted in Fig. 8.7(a), where we plot the color coded reflection 
amplitude against the probe frequency from the VNA and the applied current. We 
observe two anticrossings when all three components hybridize and obtain the 
two tripartite point at I4 = 3.909 A and In = 4.922 A with an effective coupling of 
ĝı/27 = 18.71 + 0.05 MHz and 22/21 = 23.07 + 0.08 MHz, respectively. Probing 
witha power of P = —95 dBm decouples the qubit from the system, see Fig. 8.7(b), and 
we observe the cavity magnon-polariton analog to in Sec. 7.2. The coupling strength 


of grm = 30.41 + 0.06 MHz is a factor of 2.93 larger compared to the measurements 
with the 0.5 mm diameter YIG sphere, which matches an expected 2.83 times larger 
value due to the increased number of participating spin quite well. 


We fit the tripartite resonance frequencies to the eigenvalues of the coupling matrix, 
Eq. (4.6), considering the constant resonator frequency wr to be constant and the 
qubit frequency changing linearly with the applied magnetic field, see Sec. 4.3, 
and obtain 


w,/2n (I) = 5.430 446 GHz — 171.7kHz 4 and (8.13) 


wl /27 (I) = (48742 + 0.1286 GHz. (8.14) 


The frequencies are also sketched into Fig. 8.7(a). The individual coupling strengths 
are fitted to 


841/27 = 67.49 + 0.02MHz, (8.15) 
gin /27 = 30.32 + 0.04MHz, and (8.16) 
Sqm/21t = 0.661 + 0.038 MHz. (8.17) 


The prime symbol indicates the parameters that were obtained by fitting to the 
tripartite coupling matrix, in contrast to ones obtained by coupling only two 
system components. 


The fitted values compare well to the ones obtained from the measurements 
involving only two system components. The zero-current value of the qubit 
frequency deviates from the one obtained by the power dependent spectroscopy by 
wg WG (I = 0) = 8.7kHz, which confirms the expected magnetic field dependence. 
We also recover the coupling strengths with an accuracy within an order of 100 kHz or 
< 0.5%. The fit confirms the assumption of an almost negligible direct qubit-magnon 
coupling, contributing an order of magnitude less, compared to the indirect coupling 
via the TE102 coupling mode. 
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8.3.2 Input-output fit 


We models the full reflection spectrum of the tripartite system by a fit to an analytic 
model derived using the input-output formalism, see Sec. 4.2. The data processing 
prior to the fit is analog to the procedure described in Sec. 7.2 for the magnetic 
linewidth measurements, here including the three tripartite branches. We estimate 
the field independent background signal from the measurement setup, include it 
into the fit function and reduce the number of free fit parameters to two, the qubit 
linewidth kg and the magnon linewidth km by using the previously determined 
component frequencies and resonator linewidths as constant fot parameters at the 
corresponding magnetic field values. Figure 8.8 shows the reflection amplitude 
and the corresponding fit plotted against the probe frequencies for the applied 
current values closest to the two tripartite working points at 15"** = 3.91 A (a) and 
['*8 = 4.92 A (b), together with the corresponding effective coupling strengths 21 
and 22. 
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Figure 8.8: Reflection spectrum at the tripartite working points Data and corresponding analytic 
input-output fit of the amplitude of the reflection spectrum closest to the tripartite working points at 
I = 3.91 A (a) and I, = 4.92 A (b). The fit qualitatively recovers the three main resonance frequencies, 
but the resonance line shapes are not reflected most accurately, due to a not precise enough modeling 
of the background signal. Especially the smallest resonance features are influenced strongest by any 
inaccuracies in the background estimation. This reduces also the physical meaning of the numerical 
values representing the optimal fit parameters. Nevertheless, the plots clearly show the anticrossing 
splitting of the tripartite working points and confirm the strong coupling between the three components. 


While we recover the general shape of the spectrum and the three main resonances, 
there is however a strong deviation of the fit from the data regarding the resonance 
line shape. While the most pronounced and deepest resonance are represented well 
by the fit at both currents (for visualization not completely visible in the plots), the 
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less pronounced resonances are not shaped correctly, failing to represent both the 
width and the depth of the dips correctly. We attribute this to inaccuracies in the 
modeling of the background signal, which becomes more crucial at smaller signal 
strengths. The fit fails to correctly model the only 2 — 5 dB signal from the dressed 
frequencies on top of the background. This is especially well visible at the lowest 
resonance in Fig. 8.8(a) and the highest resonance in Fig. 8.8(b). While the resonance 
in the first case is too broad and to shallow compared to the data, with the fit curve 
rising rising before the data, the fit overestimates the depth of the resonance in 
the second case, at the 5.53 GHz resonance in subfigure (b), and the linewidth is 
modeled much more narrow than in the measurement. At the same time we observe 
a shoulder feature in the fit that does not match the data, where the background is 
over estimated compared to the appearing resonance. 


Plotting the analytically obtained reflection spectrum color coded, we recover general 
course of the pronounced resonance features with increasing current, but fitting 
the data to extract linewidths from the shape of the resonances does not result in 
reliable data. In order to achieve such a precise fit, there is need for a more elaborate 
background modeling routine. 


8.3.3 Temperature dependent reflection spectra 


The weak-excitation approximation of the qubit in the input-output formalism treats 
the qubit effectively as a two-level system, ignoring all higher non-harmonic energy 
levels. While we operate the system far away from the temperature equivalent of the 
participating resonances, the temperature does not necessary reflect the electron 
temperature or the occupation number of the qubit. In order to visualize the effects 
of excited higher qubit modes in the reflection spectrum, we simulate spectrum of 
the tripartite system Hamiltonian with respect to the temperature and compare the 
features to corresponding measurements. 


We simulate the system reflection spectrum using qutip's qutip.correlation.spectrum 
function in an n = 3 dimensional Hilbert space for each system component. The 
function calculates the spectrum of the correlation function 


S(w)= f dr lim (A(t+7)B(t)) e lor (8.18) 


—o t>00 


of the full tripartite system Hamiltonian without the rotating wave approximation, 
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with the known coupling strengths and the field dependent frequencies. We construct 
the qubit Hamiltonian in the charge basis, 


2 
Ay=Ezc (2-3) mal- Fiata (820 


where we calculate the eigenenergies and multiply them with the correspond- 
ing excitation number basis. The coupling terms to the magnon and resonator 
system are calculated individually for the anharmonic qubit levels based on 
the generalized Jaynes-Cummings Hamiltonian, here stated exemplary for the 
qubit-resonator coupling, 


Agr = Lħgy |i) (jl (a +a"), (8.21) 
4] 


with the coupling energies 
gi œ (illj). (8.22) 


The proportionality factor is based on the fitted coupling strength gqr, describing 
the coupling to the first qubit transition. 


The resonator and magnon Hamiltonians and their coupling are constructed as in 
Eq. 2.38, with the corresponding creator and annihilator operators and the fitted 
coupling strength. We add dissipation to the system by defining collapse operators 
that are introduced in the Lindblad equation 


d i ry R X Er 
gô O = -g [Plays p Ol +E (Eso (NCL 3 [Een 0) » (8.23) 


where the collapse operators E, = VGA j describe the operator A j through which 
the system couples to the environment with the coupling rate xj. The interaction 
enables both the decay and the creation of excitation. We consider excitation creation 
and annihilation in the resonator and the magnon due to coupling to a bath and 
internal losses. The collapse operators of the qubit models the energy relaxation 
and thermal population at a 1/T} rate combined with an addition dephasing rate. 
The environment temperature is accounted for by attributing a thermal occupation 
number nn = exp (-hw;/kgT) based on the Boltzmann distribution to the individual 
components. It is included as an additional term in the annihilator operators and asa 
proportionality factor in the collapse operators describing an excitation. Exemplary, 
we find for the collapse operators of the photon states in the resonator 


En = yr (14 nh) a (8.24) 


Ca = y kata. (8.25) 
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Figure 8.9: Tripartite reflection spectrum and simulation (I) (b) The T = 25mK simulation shows a 
rich spectrum of thermally excited modes besides the fundamental one photon modes described by 
the analytic input-output formula. The corresponding data (a) shows the fundamental mode and an 
additional two photon mode going to higher currents. (d) The simulation at 100 mK shows a thermal 
broadening of all modes. (c) The data shows more two photon transitions that are also more pronounced 
and broadened. The first hints of three photon transitions are visible as well. The energy eigenvalues 
of all modes that are visible in the data are indicated in the simulation plots, with colors indicating 
the number of photons participating in the transition. The color gradient of the simulated correlation 
spectra is adjusted to visualize all possible transitions best. 


We measure the reflection spectrum of the tripartite system with parameters similar 
to the ones described in Sec. 8.3.1 for temperatures up to 500 mK, to study the 
influence of the thermal environment on the system and compare the data with 
the corresponding simulated spectra. In order to have a similar color, we plot the 
negative correlation spectrum to the data. This maps the high correlation to the 
resonance dips in the data. 
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Figure 8.9 shows a side-by-side comparison of the measured reflection spectrum 
on the left hand side, subplots (a) at T = 25 mK and (c) at T = 100 mK, with their 
corresponding simulation on the right hand side in subplots (b) and (d). We increase 
the visibility of the data, by highlighting the eigenenergies of the modes appearing 
in the measurements as dashed colored lines in the simulation. 
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Figure 8.10: Tripartite reflection spectrum and simulation at (II) At T = 250mK, the measured 
spectrum (a) shows already a decoupling of the qubit from the system, where several mixed modes 
are visible. We observe the magnon-polariton at T = 500 mK, without qubit influence (b). In contrast 
to these measurements, the corresponding simulations still show a qubit influence and the tripartite 
spectra is visible, but thermally smeared out, (b) and (d). 


At T = 25mkK, the simulation shows a rich mode spectrum in addition to the 
low frequency tripartite modes, which are indicated as white dashed lines in 
the simulation data. Several other magnon frequency dependent modes, which 
correspond to multi photon transitions, are visible in the simulation and we are able 
to identify one of these transitions in the measured spectrum. At high currents, we 
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find a small resonance at around 5.45 GHz, which for lower currents shifts to lower 
frequencies. This mode corresponds to a two photon excitation in the system and is 
indicated as a red dashed line in the simulation plot. Increasing the temperature 
to 100 mK, comparing the subplots in the lower row, the resonances in both the 
measured and simulated spectrum broaden and become more pronounced due to 
the increased interaction with the environment and thermal population of higher 
modes. Many of the simulated modes are visible in the measured spectrum. Besides 
the mode seen at 25 mK, other two photon modes (red dashed lines) and the onset 
of some three photon modes (blue dashed lines) are visible in the measured data 
and indicated in the simulation plot. 


Increasing the temperature further, the qubit-resonator coupling breaks down 
similar to Fig. 8.6(a) and the spectrum transforms into the known magnon-polariton 
spectrum without qubit participation at T = 500 mK, see Fig. 8.10(c), witha transition 
regime at T = 250 mK shown in Fig. 8.10(a). The qubit is not yet decoupled in the 
simulations, the smeared out modes of the tripartite spectra are still visible at higher 
temperatures, see Figs. 8.10(b) and (d). 


The qutip simulations are able to reproduce the measurements in the temperature 
regime, where the qubit is not thermally saturated quite well. We are able to reproduce 
the observed broadening of the linewidths due to the increased coupling to the 
environment and the simulations allow us to identify the observed multi-photon 
transitions at higher environment temperatures. 


8.4 Towards Rabi oscillation in the strongly coupled 
tripartite system 


Coherent energy exchange, which was shown by the strong coupling at the tripartite 
working points, between the system components is necessary for application in 
quantum information processing. In time domain measurements, this manifests 
as Rabi oscillation between the system components. Similar to the pure qubit 
characterization in Sec. 8.4, we analyze the system by applying pulsed microwave 
excitation and monitor the time dynamics to obtain excitation and coherence 
lifetimes. First, we characterize the strongly qubit-resonator system with the magnon 
frequency far detuned at zero applied current. Afterwards we ramp the magnetic 
field to the first tripartite working point and perform current dependent Rabi 
oscillation measurements. 


Rabi oscillations are measured by applying microwave pulses at various lengths 
with the previously determined qubit frequency to the system and measure the qubit 
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state. The pulse rotates the Bloch vector between the ground and excited state and 
plotting the qubit inversion against the pulse length shows the continuous change in 
qubit state with an exponential excitation decay. By fitting a damped sine function 
we obtain the length of the n pulse, Tr, that rotates the qubit state by 180°. The 
measured data together with the fit and the pulse scheme is shown in Fig. 8.11(a) 
and we obtain T = 27.14 ns for the strongly coupled qubit-resonator system. 


a b 
ir? " q F af ar 7 
TA = 27.14ns T =152.5ns 
x 
5 5 (m) A readout 
2 Y t 
wo © Pr = 
205 E WW AN 
5 5 
= = 
o o 
ob + J Ole , s 
0 100 200 300 0.0 0.5 1.0 
c Pulse length t (ns) d Delay At (us) 
5 er; rn r 7 a r p m 
T3 = 168.3 ns (3) t m% (3) readout 
c E 2 2 
5 s Mmmm AAN 
© 0.5} © 
É 205 š 
5 5 
= 3 
o o readout| O 
; MWE NW i TI =221.3ns 
¿E , i i OL ee , i i , 4 
0 100 200 300 0 100 200 300 400 500 
At(ns) Delay At (ns) 


Figure 8.11: Dissipative dynamics of the qubit-resonator mode (a) Driving the qubit transition 
frequency for varying pulse lengths t rotates the qubit Bloch vector on the sphere between the ground 
and excited state. Exciting the qubit takes a time T? = 27.14ns. (b) Measuring the qubit state after 
excitation and subsequent time delays results in an energy relaxation time T = 152.5ns. Decoherence 
measurements in the Ramsey sequence (b) and with an additional Hahn echo pulse (c) result in qubit 
coherence lifetimes T;’" = 168.3ns and TF = 221.3ns, respectively. No current was applied during 
these measurements. 


We measure an energy relaxation time TE = 152.5ns and a decoherence time by a 
Ramsey pulse sequence Te = 168.3 ns, see Figs. 8.11(b) and (c). The corresponding 
pulse sequence is sketched in the corresponding plots and described in more detail 
in Sec. 8.4. Additionally we perform a Ramsey type measurement with an additional 
Hahn echo pulse [201], a n pulse that refocuses the dephased qubit state, after 
half the delay time. We measure an spin echo decoherence time T} = 221.3 ns, see 
Fig. 8.11(d). All measured lifetimes are listed in table 8.3. 
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8 Magnon based nonlinear quantum system 


Comparing theses lifetimes to the transmon lifetimes in the dispersive limit we find 
the time of the n pulse a factor 5.9 smaller in the strongly coupled case. The smaller 
detuning between qubit and resonator allows for a faster energy transfer from the 
microwave feed line to the qubit through the resonator bandpass. The T; time and 
the T% time are smaller by a factor of 2.6 and 3.1, respectively, which we attribute to 
an increased loss out of the resonator. 


In order to observe coherent excitation exchange between in the tripartite system, 
we sweep the current up to the first tripartite working point and perform Rabi 
oscillation measurements, shown in Fig. 8.12. Figure 8.12(a) shows the color coded 
qubit inversion with respect to the applied pulse length and the coil current. The 
individual Rabi oscillation measurements at the two current values marked by white 
dashed lines at 3.45 A and 4.07 A, with the latter corresponding to the current value 
of the first tripartite point are show in subplot (b). The current value of the first 
tripartite point is slightly different compared to the one in Sec. 8.3.3, since these 
experiments were performed in a separate measurement run. 


The qubit inversion shows a decreasing pulse length for the first qubit inversion 
maximum with increasing current but the overall signal strength decreases, as seen 
by the data colored all in blue. While we still measure Rabi oscillations with a 
maximum qubit inversion of 75.6 % at 3.45 A, the signal strength decreases with 
higher current and the Rabi oscillation vanishes completely at the tripartite point. 


Table 8.3: Dissipation parameters of the strongly coupled qubit-resonator system The table includes 
the duration of the n pulse T¥, the energy relaxation time T}", and the decoherence times obtained by 
a Ramsey measurement with and without an additional Hahn echo pulse T and T) "Y, respectively. 
No current was applied to the coil at these measurements. 


Tr (ns) TE (ns) TF (ns) Te (ns) 


27.14 152.5 168.3 221.3 


We are able to implement nonlinear hybrid systems using magnon excitation. In 
the dispersive regime we observe indirect qubit-magnon coupling and observe the 
square root dependence of the coupling strength on the number of participating 
spins. While the effective coupling places the system in the strong coupling regime, 
the signal to noise ratio is too low to observe the hybridized mode dips in the 
spectrum. We increase the coupling by reducing the detuning in the strongly coupled 
tripartite system, where we observe and fit the three dressed modes. Thermal 
excitation of higher system transition are observed and simulated as well. Although 
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8.4 Towards Rabi oscillation in the strongly coupled tripartite system 


the system shows strong coupling in the spectroscopic measurements, we do not 
observe Rabi oscillation in the system approaching the first tripartite working point. 
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Figure 8.12: Current dependent Rabi oscillations (a) Color coded qubit inversion obtained by Rabi 
oscillation measurements with respect to the applied coil current. The maximum qubit inversion 
shifts to lower pulse lengths with increasing coil current, but the overall oscillation signal decreases 
and vanishes approaching the first tripartite working point. The color gradient is adjusted for better 
visibility. (b) Individual Rabi oscillation measurements for I = 3.45 A (blue) and I = 4.07 A (green), 
that highlight the decreasing Rabi oscillation signal. Both currents are also indicated as dashed lines in 
the color plot. The applied currents do not coincide with the current axes in Sec. 8.3 since the time 
domain measurements were measured in a different measurement cool down with a slightly changed 
current-to-field relationship. 
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9 Conclusion 


The goal of this thesis was to investigate magnetic excitation in quantum circuits, 
with the goal to explore the potential of magnon based technology for applications in 
quantum computing and quantum information processing. We embedded magnons 
in a linear and nonlinear quantum hybrid system in different coupling regimes 
and studied the microwave properties of the magnetic resonance at millikelvin 
temperatures and single photon excitation powers. The presented experiments 
demonstrate the dominant loss mechanism of magnon excitation at the working 
environment for superconducting quantum circuits and the first measurements on a 
strongly coupled light-matter-spin tripartite system. 


In the first part of the this thesis we designed a harmonic photon-magnon hybrid 
system still in the regime of classical physics. We strongly coupled the Kittel mode 
magnons in a ferrimagnetic YIG sphere to cavity resonance photons, creating a 
cavity magnon polarion quasiparticle with a coupling strength g = 10.39 MHz. We 
extracted the internal magnon resonance linewidth from the microwave reflection 
spectrum and identified the limiting loss mechanism in different temperature 
regimes, with special focus on the millikelvin temperature regime in a dilution 
refrigerator. Together with our collaborators at JGU Mainz, we confirmed scattering 
on rare-earth impurities in the YIG and multi-magnon scattering on the sample 
surface as the dominant loss mechanism when cooling the system from room 
temperature down to single digit kelvin temperatures [76]. Cooling further to 
millikelvin temperatures, the magnetic resonance linewidth is limited by excitation 
losses into two-level system, a characteristic loss model observed in glasses and 
superconducting circuits as well. We confirmed the temperature trend of the 
linewidth following a hyperbolic tangent and observed the power dependence of 
the linewidth following a 1/ vP trend from single photon to the 107 exchtation 
regime, which is attributed to incoherent losses to a bath of two-level systems as 
the main loss mchanism. With our measurements we expanded the study of cavity 
magnon polaritons to off-resonant measurements, where we were able to extract the 
magnon excitation percentages in the cavity magnon polariton. With the number 
of excited magnons depending on the frequency detuning between magnon and 
photon at constant experimental input power, the observed magnon linewidth trend 
resembled the power dependent measurements [62]. 
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9 Conclusion 


Having demonstrated the strong magnon-cavity coupling and usability at low 
temperatures, we fabricated a superconducting 3D transmon qubit as a nonlinear 
element for our quantum circuit. This allowed us to address single magnon excita- 
tions, in contrast to the Boltzmann-distributed excitations of the magnon-photon 
system. We demonstrated the quantum properties of the qubit by Rabi oscillations 
and measured quantum excitation lifetimes on the order of 500 ns. With the system 
in the dispersive limit, the magnon frequency is tuned towards the qubit and 
magnon-vacuum induced Rabi oscillations as an anticrossing in dispersive shift 
measurements were observed. The Qubit and the magnon coupled via virtual photon 
exchange of the surrounding cavity modes with an effective coupling strength of 
3.3 MHz, which did not put the system in the strong coupling regime. Measure- 
ments on an improved system with a larger coupling of 9.34 MHz due to a larger 
number of participating spins showed a signal to noise ratio not strong enough for 
further investigations. 


As a next step, we built a quantum hybrid system that coupled light, matter, and spin 
excitation to increase the effective coupling strength and eliminated the need for 
dispersive shift measurements. This allowed us to measure the system directly via 
microwave spectroscopy. We aimed at the cavity resonance modes to be degenerate 
with the qubit transition, reaching a detuning of 25 MHz compared to a 67.6 MHz 
coupling strength. Tuning the Kittel mode frequency through the already dressed 
states revealed two tripartite working points with strong effective coupling on the 
order of 20 MHz between the components. We were able to fit the reflection spectrum 
with the input-output formalism and simulated the spectrum of the full system 
Hamiltonian, which allowed us to identify multi-photon transitions. Temperature 
sweeps and corresponding simulation showed thermal excitation and subsequent 
transitions containing up to three photons, before the qubit functionality broke 
down due to thermal effects. 


In the final experimental section we investigated the dissipative dynamics of the 
tripartite system, looking for Rabi oscillation between the system components. 
Without magnon contribution, we observed relaxation and coherence times about a 
factor of 3 smaller than for the bare qubit in the dispersive regime. Increasing the 
Kittel mode frequency however, the Rabi oscillation frequency decreased slightly 
before the signal vanished completely towards the tripartite point. 


In conclusion, we identified the dominant magnon loss mechanism at low temper- 
atures and presented coherent exchange between a superconducting qubit and a 
Kittel mode magnon in the dispersive regime, as well as proof-of-principle mea- 
surements on a strongly coupled tripartite system, a nonlinear quantum circuit 
containing magnons. 
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9 Conclusion 


The presented experiments are a promising foundation for more sophisticated 
measurements, for which the system could benefit from modification or redesign of 
individual components. 


Outlook 


During the course of this thesis, a study showing the first application of a quantum 
magnonic system besides a proof-of-concept experiment was published, where the 
number of excited magnons was expicitly counted [67]. This experiment is based 
on coupled qubit-photon systems measured more than 10 years ago [169, 170] 
and shows one roadmap for quantum magnonics, realizing features previously 
measured with photons using magnon excitation, with possible next candidates 
being the magnon equivalent of coherent superposition states [171] and Schrédinger 
cat states [172]. A recent preprint follows this path in creating the single-photon 
detector counterpart for a magnonic system [68]. 


Time resolved measurements regarding Rabi oscillation or direct measurements of 
magnon lifetimes, besides an estimation based on a resonance linewidth are yet to 
be shown. This would open up the door towards quantum storage experiments, 
where a quantum state is swapped to the magnon and back to the qubit to be read 
out again. Quantum gates via parametric coupling or fast qubit detuning via as 
Stark shifts are already available in the toolbox for superconducting circuits and the 
measurement protocol can be shortened by single shot qubit readout with active 
reset. Experimentally further away are quantum transport measurement, where a 
quantum excitation is encoded in k 4 0 magnons that propagate through the sample 
to be picked up by a second qubit. 


Looking at our current measurement setup regarding the next experimental gener- 
ation, there is room for improvement on all three system components. A manual 
polishing treatment of the YIG surface may reduce the impact of TLS losses and 
create a sharper magnon resonance linewidth with increased excitation lifetimes. 
Conference presentation and comments by members of the Tang lab at Yale University 
suggest a much improved resonance linewidth after a manual ultra-fine polishing 
of the sample surface, compared to the values obtained with the initial surface 
treatment by the supplier. The cavity resonator design can be improved regarding 
coupling strength of to the Kittel mode, where a coupling in the range of 50 MHz 
would be advantageous. Fabricating the cavity out of aluminium instead of copper 
may improve the Q factors, despite the superconductivity breaking down in some 
parts due to the applied magnetic field. At last, an improved qubit fabrication leads 
to longer coherence times which increase the overall experimental possibilities. 
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9 Conclusion 


With these possibilities and already available approaches to solve similar problems 
we are convinced that magnons can play a role as a building block in future 


quantum circuits. 
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Appendix 
A Fabrication parameters 


Table 1: Substrate cleaning 


solvent ultrasonic bath IPA HO hot plate 
NEP 10 min yes yes 100°C,5min 


Table 2: Application of the resist stack 


resist pre spin reflow hotplate 
LOR 500nm 300rpm, 10s 3000 rpm, 60s 30s 180°C, 5min 
PMMA 950K 


AR-P672,03 > Orpm,10s 6000rpm,60s 30s 145°C,5min 


Table 3: Electron beam lithography with a 10 keV JEOL writer 


writing current dose 
STRUCTS 1nA 300 pC cm? 
STRUCTS 100 pA 300 pC cm~? 
STRUCTS 100 pA 60 pC cm? 


115 


Appendix 


Table 4: Development 


developer time rinse dry on spin coater 
MIBK/IPA 1:3,25°C 65s IPA (655), H2O (60s) 7500rpm, 120s 
MIF 726/H20 3:2 120s H20 (60s) 7500 rpm, 120s 


Table 5: Plasma clean with Kaufman source, PLASSYS 


as mixture lasma voltage current time 
8 P 8 


Ar (4sccm), O2 (0.5 sccm) 200 V 10mA 20s 


Table 6: Al evaporation, shadow angle, PLASSYS 


angle rate thickness 
1 


Layer 1 8° 0.2nms” 30nm 


Layer 2 8° 0.2nm s7! 50 nm 


Table 7: Lift-off 


solvent temperature time ultrasonic bath rinse 


NEP 90°C 4h 1min IPA, H20 
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quantum systems of different physical origin, with electromagnetic excitation 
being used to connect the individual components. While the underlying logic 
with fast and scalable quantum gates is often realized using a solid-state ap- 
proach utilizing i.e. superconducting circuits, additional elements are needed for 
storing quantum information or frequency transducing. These circuits operate 
at microwave frequencies and allow for research of fundamental interactions 
of physical systems with light. With the nonlinear energy structure of a super- 
conducting quantum bit, the interaction of matter with other physical systems, 
i.e. spins, can be studied in similar fashion. In order to achieve strong coherent 
coupling of superconducting circuit elements to magnetic spin excitation the 
number of participating spins is increased and collective magnetic bulk excita- 
tion are employed to study spin-based hybrid systems. 


In this book, hybrid systems based on yttrium-iron-garnet (YIG), three dimen- 
sional microwave cavity resonators, and superconducting transmon qubits, are 
investigated by continuous wave and pulsed microwave spectroscopy. Limitations 
to the magnetic linewidth in the quantum regime are identified and coherent 
exchange between a magnon and a superconducting qubit are demonstrated. 
Finally, a first step towards a strongly coupled hybrid system containing all three 
components is demonstrated. 
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